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ABSTRACT 

This  report  covers  the  efforts  undertaken  to  develop  a 
mathematical  model  of  the  FAMMS  I  surf ace-to-surface  missile 
suitable  for  preliminary  analysis  of  the  effects  of  digital 
computer  sampling  on  system  performance.  Planar  equations  are 
developed  and  linearized  to  permit  digital  autopilot  design 
to  procede.  A  design  methodology  is  developed  and  described 
through  application  to  a  simplified  version  of  the  FAMMS  I 
equations. 
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SECTION  I.  INTRODUCTION 

The  purpose  of  this  work  is  to  establish  the  technical  ground¬ 
work  for  the  analysis  of  the  performance  of  a  digital  guidance  and 
control  (G&C)  system  for  use  on  an  advanced  Future  Army  Modular 
Missile  System  (FAMMS).  Presently  FAMMS  is  envisioned  to  be  a 
surface-to-ai r  missile.  The  tasks  described  to  achieve  this  analysis 
are  envisioned  as: 

1.  Derive  equations  of  motion; 

2.  Linearize  equations  of  motion; 

3.  "Close  the  loop"  by  including  guidance  and  control 
(autopilot); 

4.  Design  and  stabilize  autopilot; 

5.  Develop  simplified  planar  time-varying  simulation; 

6.  Develop  candidate  target  models; 

7.  Design  guidance  loop(s); 

8.  Conduct  performance  studies. 

The  work  reported  upon  in  this  document  comprises  Tasks  1-4 

as  described  below: 

1.  Derive  planar  equations  of  motion; 

2.  Linearize  planar  equations  of  motion; 

3.  Close  the  loop  with  a  digital  autoDilot,  assuming 

a  "perfect"  guidance  command;  j 

4.  Design  and  stabilize  a  postulated  autopilot. 

The  performance  of  the  postulated  autopilot  is  checked  by 
application  of  the  "Cross-Multiplication"  Technique  described  in 
Appendix  A.  It  is  anticipated  that  the  simulation  philosophies 


2 


described  in  Appendix  B  will  be  used  in  part  or  as  a  whole 

for  detailed  performance  studies  with  the  MICOM  FAMMS  Simulation 
presently  under  development. 


j 

t 
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SECTION  II.  EQUATIONS  OF  MOTION 

The  planar  equations  of  motion  of  the  FAMMS  were  derived 
considering  the  three  coordinate  systems  as  depicted  in  Figure  1 . 
The  three  coordinate  systems  are  defined  as  folTows: 

A  A 

1R,  3R  -  U  it  vectors  defining  an  inertial 

coordinate  frame  used  for  reference, 

A  A 

lb,  3b  -  Unit  vectors  defining  a  body  fixed 
(missile)  coordinate  frame, 

A  /\ 

IT,  3T  -  Unit  vectors  defining  a  target  fixed 
coordinate  frame. 

The  lb,  3b  rotates  through  an  angle  6  with  respect  to  the  1R,  3R 
inertial  reference  system.  The  origin  of  this  body  fixed  system 
is  defined  by  the  vector  Rm  ,  the  engine  gimbal  angle  is  defined 
by  the  angle  6^,  and  the  aerodynamic  surface  deflection  by  the 
angle  5.  The  vector  equations  of  motion  are  then  as  follows: 

‘ '  A  A 

MR^  =  (F-A)  lb  +  ( F  5  - N )  3b  +  mg  3R  ( 1 1 - 1 ) 

10  «  M  -  t0N  -  l  FSe  ( 1 1 - 2 ) 

Where 

F  =  Engine  thrust  force 
A  =  Aerodynamic  axial  force 
N  =  Aerodynamic  normal  force 
M  =  Aerodynamic  pitchinq  moment 

l0  =  Distance  from  vehicle  center  of  mass  to  moment 
reference  point 

l  =  Distance  from  vehicle  center  of  mass  to 
^  engine  gimbal  point. 
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Now  representing  the  vector  as  follows: 

R  =  X  1R  +  Z  3R, 

— m 

yields  the  following  scalar  equations  in  the  inertial  coordinate 
frame 


mX  =  (F-A)  cos0  -  (F6  -N)  sin9 

( I 1-3) 

mZ  =  (F-A)  sin0  +  (F5  -N)  cos9  +  mg. 

( 1 1-4) 

Equati ons  ( 2 )  ,  (3),  and  (4)  represent  the  three  nonlinear  and  time 
varying  differential  equations  that  must  be  solved  to  yield  the 
vehicle  motion  as  a  function  of  time.  The  autopilot  and  guidance 
system  will  provide  <5  and  <5^  as  a  function  of  the  error  conditions 
The  equations  can  then  be  numerically  integrated. 

For  stability  analyses,  a  linear  set  of  equations  is  required 

/V  A 

defined  in  a  body  fixed  reference  system  ( 1 b -  3  b ) .  To  derive  these 
linear  equations  the  following  definitions  are  made 
r  =  1  9  N  I 

Na  -  JIT  3a  !a  =  0  ( 1 1-5) 

r  _  1  3Mn  ! 

m  “  SQD  3a  I  (II-6) 

'  a  =  0 


C 

Aa 


1 

~S$ 


3A 


3a 


a  =  0 


( 1 1-7 ) 


C 


N6 


1  3N 

SQ  96  ■ 6=0 


( 1 1-8) 


6 


^m6 

1 

SQD 

3M 

36 

■  ( 

5=0 

(I 1-9) 

CA6  S 

1 

"Try 

3A 

36 

a=0 

(11-10) 

^md  = 

i 

3  M 

SQD 

3  a 

a  =  0 

(11-11) 

where 

D  = 

Reference  vehicle  diameter 

S  = 

ttD^/4 

Q  = 

(H)pV2 

Thus  the  linear  aerodynamic  forces  become: 


N  = 

55  <CN0“ 

+  CN66)  +  N0 

(11-12) 

M  * 

DSQ(Ctoc 

+  Cm56  +  C»4&)  +  M0 

(H-13) 

A  = 

SQ  (CAa« 

+  CA55)  +  A0 

(11-14) 

The  left  hand  side  of  Equation  1 1  - 1  can  be  transformed  to 
a  body  fixed  system  through  consideration  of  the  following 
definition 

=  Ulb  +  W3b  (11-15) 

=  lb  component  of  velocity 

A 

=  3b  component  of  velocity. 

V  =  (U  +  0W)Tb  +  (W  - 8 U )  3b  (11-16) 

=  UQ  +  u(t) 

=  w(t) 

*  s0  +  e(t) 


where 


Thus 


Now  let 


u 

w 


9 


(11-17) 


The  quantities  u,  w  and  0  represent  small  deviations  from  the 


nominal  flight  condition  represented  by  UQ  and  Qq. 
noted  that, 

w 

“  *  h  ' 

Further,  define 


It  is  also 


(11-18) 


(11-19) 


The  three  governing  equations  of  motion  can  be  written  as  follows 
considering  the  vehicle  to  be  initially  in  static  equilibrium. 


Mu' 

.  QS 
'  uo 

(cAa“  *  CajJ>  *  §°S  °'  9 

(11-20) 

M  a 

it 

3 

CD  • 

-  <CNca  +  CN{{)-  DJ  9 

(11-21) 

I  0 

=  oso 

•  ^ CMaa  +  CM6<5^  "  £0QS^CNaa  +  CN56^ 

+ 

OSD  C_ »a 
ma 

(11-22) 

The  further  assumption  is  made  in  these  equations  in  that  the 
only  control  variable  is  the  aerodynamic  fin  deflection  angle, 5 . 
The  design  of  the  autopilot  can  be  facilitated  in  this  study 
through  consideration  of  Eqs.  (11-21)  and  (11-22)  only  since  there 
is  typically  weak  coupling  with  the  translational  Eq.  (11-20). 
(In  aircraft  terminology  this  is  the  ’’short  period"  mode.)  If 
it  is  further  assumed  that  the  missile  is  in  level  flight  the 
transfer  function  between  vehicle  rate,  0,  and  fin  deflection, 

<5  i  s 


9(S) 
5  ( s  ) 


(s2  +  cNcs  -  CMa) 


(11-23) 


where 


(11-24) 


N6 

5 

N6 

(11-25) 

CM<5 

i  'CMS 

-  —  C  ) 
0  '"NS 

DOS 

I 

(11-26) 

CMa 

5  (CHci 

-  --0-  C  ) 
D  Na 

DOS 

I 

(II-27) 

Now  for  a  design  example  take  the  following  flight  conditions 


Mach  number  =  5.0 

Altitude  =  10,000  meter 

0O  =0  degrees 

The  FAMMS  vehicle  characteristics  at  these  flight  conditions 
are  as  follows 


CMa  •  -1-7965 

CNa  «  .064 

CM6  =  -.381 


C 


N<5 


.002 


The  required  values  for  the  transfer  function  (11-23)  are  therefore 


CM<5 

^Na 


-15.9 

.0006 

.0002 


-65.92 
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SECTION  III.  DIGITAL  AUTOPILOT  DESIGN 

It  is  assumed  that  the  task  of  the  digital  autopilot  is  to 
match  the  actual  angular  rate  of  chanqe  of  the  missile  to  the 
commanded  rate-of-change .  In  Section  II  these  variables  are 
denoted  9  and  9C,  respectively.  The  autopilot  attempts  to  achieve 
this  goal  by  developing  a  commanded  aerodynamic  vane  deflection 
angle,  S,  which  is  a  function  of  the  error  between  the  command 
angular  rate-of-change,  $c,  and  the  sensed  angular  rate-of-change, 

9.  It  is  assumed  that  the  sensed  variable  is  exactly  equal  to  the 
actual  variable  and  that  the  commanded  aerodynamic  deflection  is 
exactly  equal  to  the  actual  deflection.  These  assumptions  can 
be  relaxed  in  subsequent  more  detailed  design  efforts. 

In  order  to  cope  with  expected  forms  of  0C,  a  simple  digital 
fliter  form  is  used  for  the  autopilot  controller,  D(z).  The 
system  is  described  in  block  diagram  form  on  Fig. 2  .  The  fol¬ 

lowing  definitions  are  used: 

e  «  9C  -9,  ( 1 1 1-  1 ) 

0  1  -  e'sT  ("Zero  Order  Hold"),  ( 1 1 1 - 2 ) 

bh0  s  "  s 

and 

G  (s)  _  9  ( s  )  _  hs _  (developed  in 

P  9  O 

5(s)  s^  +  u>  Section  II).  (HI-3) 


Development  of  Autopilot  Equations. 


The  autopilot  equations  are  developed  by  applying  SAM 


(Seltzer's  Simplified  Analytical  Method)  to  the  system  described 
in  Fig.  l.*  The  six  states  (X.,  X^,  ....  X^)  may  be  assigned  by 
the  designer.  They  are  shown  in  Fig.  2  and  defined  below  in 


Column  1. 


n 


TABLE  I.  SAM  EQUATIONS 


Oriqi nal 

Eauations 

Modified  Equations 

Pul sed 

1  Equations 

X1  = 

• 

e  = 

•  • 

ec-e 

V 

S 

ec*  -  x6* 

x2  = 

xi* 

X2* 

s 

X  * 

*1 

II 

CO 

X 

D  X2 

X3  *  °  V 

X  * 
*3 

= 

D  *  Xj* 

x.  = 

X  * 

X  * 

— 

X,* 

4 

3 

4 

3 

X5  = 

Gh0 

X4 

**  *  "ho  X3* 

X  * 

a 

So*  v 

X6  = 

9  = 

GP  X5 

S'SSo  V 

X  * 
6 

a 

9*  =  GPGhO 

Applying  the  rules  developed  under  SAM,  state  equations  X-j,  X^, 

and  X,  must  be  modified.  These  modifications  are  shown  in  the 
b 

second  column  of  Table  I.  The  pulsed  equations  are  shown  in  the 
third  column.  As  described  in  Ref.  1,  the  asterisk  is  used  to 
symbolize  the  samping  operation,  be  it  real  or  analytical.  The 
pulsed  concept  may  be  simply  described  as  providing  the  value  of  the 
pulsed  variable  at  each  samping  instant  T. 

B .  Development  of  Closed-Loop  Transfer  Function. 

The  pulsed  equations  developed  above  may  be  used  to  determine 
the  values  of  the  various  system  variables  at  (or,  with  suitable 
modifications  described  in  Ref.  1,  between)  sampling  instants.  How¬ 
ever,  no  indication  is  easily  obtained  of  system  stability  or  desired 
performance.  One  way  of  obtaining  these  two  important  measures  is 

to  develop  the  closed-loop  transfer  function  of  the  system.  This 

»  • 

is  done  by  solving  X  *  (which  represents  8*)  in  terms  of  9  : 

6  c 


x  * 
*6 


9*  =  V 


Sp^hO* 


D*  X 


* 

1 
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<6*  *  *FhO*  D*  (V  “  x6*} 


6PGhO*  D* 

1  +  °* 


0  * 
c 


( 1 1 1-4 ) 


Equation  ( 1 1 1 - 4 )  may  be  readily  transformed  into  the  z-domain  to 
yield  the  closed-loop  transfer  function. 


e(z)  _  GpGhp(z) 


0c(z) 


1  +  GpGh0(z)  D(z ) 


( II 1-5) 


sT 


The  common  definition  of  z  is  used,  i.e.  z  =  e  ,  with  s  being 


defined  as  the  Laplace  operator,  s  =  a  +  i  co  =  -z,  +  /  1  -  rj- . 

ojn  n 

Making  use  of  Eqs.  ( 1 1 1  -  2 )  and  (III  —  3)  •  one  may  derive  the 


z  form  GpGho  (z): 


SPGhO  (2» 


a 


,  -  sT 


u2)  (s) 


i 


(bsinuiT)  (z-1) 
to  (z2  -  2z  cos  w  T  +  1) 


( 1 1 1 -6 ) 


The  form  of  digital  filter  is  the  topic  of  additional  study 
to  be  oerformed  at  a  later  date.  However,  the  design  philosoDhy 
may  be  developed  by  assuming  a  ratio  of  polynomials  (in  z-domain) 
for  a  simple  filter  that  represents  the  rotational  autopilot.  That 


simple  form  is  postulated  to  be, 

,2 


D(z)  = 


K2  +  K1  2  +  K0 
(z  -  l)2 


( 1 1 1-7) 
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Combining  Eqs.  ( 1 1 1 - 6 )  and  (III-7),  one  may  obtain 

(b  sin  wT)  (K-z2  +  K.z  +  KJ  ( 1 1 1-8 ) 

<=.6  <z)  D(z)  *  - - 2 - - - a -  • 

co  (zd  -  2z  cos  coT  +  1)  (z  -  1) 


One  may  now  substitute  Eq.  ( 1 1 1 - 8 )  into  Eq.  ( 1 1 1 - 5 )  to  obtain 
the  system  closed-loop  transfer  function, 


9(z )  _  (  b  sin  coT)  (K?z^  +  Kiz  +  Kn) 

Sc<Z)  J.0(dJ  K°  +  fJKi  *  V 

where  the  coefficients  of  the  denominator  are  defined  in  Table 


( II 1-9) 


TABLE  2.  CHARACTERISTIC  EQUATION  COEFFICIENTS 


j 

d  . 

J 

f . 

J 

9j 

0 

b  sin  coT 

0 

-0) 

1 

0 

b 

sin  wT 

a)  ( 1  +  2  cos  coT) 

2 

0 

0 

<2b  sin  a)T  -  to ( 1  +  2 

cos  wT) 

3 

0 

0 

0) 

The 

denominator  of 

Eq. 

.  ( 1 1 1  -9 ) , 

when  set  equal  to  zero, 

consti tutes 

the 

characteri sti c 

equation  of 

the  system. 

C. 

Selection 

of 

Numeri cal 

Values  For  Autopilot  Coefficients 

The  numerical  values  of  the  autopilot  coefficients  are 


selected  with  the  use  of  (of  course!)  the  Parameter  Plane  technique. 
The  format  of  Table  2  was  selected  to  facilitate  that  use.  The 
numerical  values  determined  in  Section  II  are  used  for  the  values 


of  b  and  u>  (characterizing  the  operating  point  about  which  the  open- 

loop  system  dynamics  were  developed).  For  this  example,  T  is 

selected  as  0.1  sec.  and  is  selected  as  0.01.  The  Control 

Dynamics  program  for  using  the  Parameter  Space  Technique  was  used. 

The  output  is  shown  in  Table  3  (following  pages).  The  A,  B,  and 

F  matrices  shown  at  the  top  of  Table  3  are  the  transposes  of  the 

d.: ,  f.,  and  g  columns  of  Table  2  when  the  above-mentioned  numerical 

J  J  j 

values  for  b,u>,K  ,  and  T  are  used. 

First,  the  coordinates  for  the  parameter  plane  are  selected 
to  be  Kg  (designated  as  PI  in  the  computer  printout)  and  (P2). 

The  stability  boundary  associated  with  the  singularity  at  z  =  +1 
is  found  to  be  (see  printout) 

P2  *  Kj  *  -  PI  -  0.01  =  -K  -  0.01  .  ( 1 1 1 - 1 0 ) 

Similarly,  the  boundary  associated  with  the  singularity  at  z  =  -1 
is  found  to  be 

P2  =  Kj  *  PI  +  4.7599  =  KQ  +  4.7599.  (  1 1  1  - 1  1  ) 

The  boundary  associated  with  the  unit  circle,  z  =  e1wnT,  is  printed 

in  Table  3  as  a  function  of  the  argument,  wnT  .  Opposite  each 

incremented  value  of  oi  T  is  Drinted  the  Jacobian  and  the  associated 

n 

values  of  PI  and  P2.  Since  the  Jacobian  (J)  does  not  change  sign, 
there  is  no  value  of  o)nT  where  a  boundary  associated  with  J=0 
occurs.  The  three  boundaries  are  sketched  on  Fig.  3  and  plotted 
on  Fig.  4  .  Finally,  a  contour  of  s  =  0.5  is  determined  and 

plotted  on  Figures  4  and  5.  The  scale  of  the  region  of  stability 


TABLE  3.  COMPUTER  PRINT-OUT  FOR  PARAMETER  PLANE  ANALYSIS 


Z  TRANSFORM  METHOD 
THE  SYSTEM  ORDER  IS  =  3 


ELEMENTS  OF  A  MATRIX  A < 0 > , A < 1 > , A < 2 > , . . . 
-1.1542E+01  O.0O00E+O0  0.00OOE+0O  0.000OE+00 


ELEMENTS  OF  B  MATRIX  B<0) , B< l > , B<2> , . . . 
O.0000E+O0  -1.1 542E+0 1  O.0000E+00  0.0000E+O0 


ELEMENTS  OF  F  MATRIX  F <  O ;■ ,  F <1 , ,  F <.  2 >  ,  .  .  . 

-8. 1 139E  +  80  i . 9293E+01  -  1 . 9403E+O1  3.  1 139E+O0 


ZETA  =  0.0O00E+00 

THE  STABILITY  BOUNDARIES  FOR  Z*+l  AND  -1  ARE  AS  FOLLOWS 
FOR  Z=-l 

P2  =  1.0000E+0O  *P 1  +  4. 7599E+00 

THE  AXES  INTERCEPTS  ARE  AS  FOLLOWS: 

AT  Pl=6  ,  P2=  4 . 7599E+00 
AND  AT  P2=0  ,  P1--4. 7599E+00 


FOR  Z*+  1 

P 2 * - 1 . 8O00E  +  00  *P1  +  -1.0000E-02 

THE  AXES  INTERCEPTS  ARE  AS  FOLLOWS: 
AT  P1*0  ,  P2»- 1  - 0000E-02 
AND  AT  P2*0  ,  Pl=-1 . 0000E-02 


INITIAL  OMEGA  =  0.00O0E+00 
NUMBER  OF  STEPS*  50 
W  n  T  =  1 . 000OE-25  DET  * 

WriT *  ©.2332E-32  DET® 

WnT*  l . 256SE-0 1  DET= 

WriT  *  1 .  SS50E-0 1  DET 

WnT*  2. 51 3SE-01  DET 

W  r*i  T  =  3.  141SE-01  DET 

WnT  =  3 • 7S99E-0 1  DET 

WnT®  4 . 3932E-0 1  DET 

WnT*  5 . 0265E-0 1  DET 

WriT»  5.S54SE-01  DET 

WnT-  b . 2332E-0 1  DET 

WnT*  S . 91 1 5 E - 0 1  DET 


FINAL  OMEGA 

IE -23 

3 . 3  b  4  5  E  ♦  0  0 
1  .  So9bE  *0  1 

2 . 4  9  6  2  E  *  0 1 
3. 3129E-01 
4 . l 1S5E+01 
4  .  ?  O  3  9  E  +  0  1 
5 .  b  7  2  0  E  +  ©  1 
S. 41 76E+0I 
" . 1 379E+0 1 

7 . 5  3  0 1 E  +  0  1 
9. 4?14E-0I 


141SE+00 

4 . 1  3  7  b  E  -  0  1 
4. 1 430E-0 1 
4,01 0SE-0 1 
3 . 7953E-0 1 
3.  5  O  4  0  E  -  0  1 
3. 14F7E-01 
'362E-01 
■374E-01 
1  7  2  E  -  0  1 
445E-01 
9  3  0  E  -  0  2 
29EE-02 


P  1  *  - 4 , 2 3 7 6 E - 0  1  Pi 

Pl«-4. 2598E-01  P: 

°l=-4. 1767E-01  Pi 

P  1  —  —  4 »  0  3  y  4  E  —  u  1  P 

*1  =  -  3 . 345SE-0 1  P- 

1  ®-  3.  :  990E-O1  F  ■ 

F 1 *-  ! . 299EE -O 1  Pi 

P  1  *  -  2  .  ?  4  3  b  E  -  0  1  ci 

P 1 = - 2 . 5  4  7  4  E - 0 1  ° 2=  1 

P l  =  - 2 , 09  7  5E - 0 1  P2=  1 

P 1 =- 1 . O007E-O1  P2  =  3 

P2--1.0590E-01  F  2=  4 


WnT*  7.5393E-01 
WnT *  3.1881E-01 
WnT*  3.79S5E-01 
WnT»  9.4248E-01 
WnT*  1.0053E+00 
WnT  *  1 . 063 1 £  +  00 
WnT  =  1.1310E+00 
WnT*  1.  1938E+00 
WnT*  t  •  2566E+00 
WnT -  1 . 3 1 95E+00 
WnT*  1.3323E+00 
WnT  *  1.4451E+00 
WnT*  1 • 5030E+90 
WnT*  1.5708E+00 
WnT*  1.6336E+00 
Wn  T  =  1  * 69o5E+00 
WnT  =  1.7593E+00 
WnT-  1 , @22 1 E+00 
WnT  =  1.8350E+00 
WnT®  t . 9478E+00 
WnT*  2.0106E+0© 
WnT*  2*  0733E+0Q 
WnT  =  2 . 1 3  o  3  E  +  0  0 
WnT®  2.  1991  £-#-00 
WnT =  2 . 2  6 1  9  E  +  0  0 
WnT  =  2.3248E+8Q 
WnT =  2.387*£+00 
WnT*  £ , 4504E+0O 
WnT*  2*51 33E+00 
WnT*  2.57SIE+00 
WnT*  2 . 63S9E+00 
WnT*  2.701 3E+00 
WnT*  2.7b46E+00 
WnT®  2 . 8274E+ 0 0 
WnT®  2 . 89036+00 
WnT*  2.953 1 £+00 
WnT*  3.01 59E+00 
WnT*  3 , 0733E+00 
WnT*  3. 1418E+00 


BET* 

9. 1 1916+01 

BET* 

9 . 71086+01 

SET* 

1 . 8264E  +  02 

BET* 

1 . 0777E+32 

BET* 

l .  12436+32 

BET* 

1  .  16746+02 

BET* 

1 . 23546+02 

BET* 

1 . 2336E+3*! 

BET* 

1 . 26696+3 2 

BET* 

1 . 2903E+32 

BET* 

1 . 3385E+02 

BET* 

1 . 3216E+32 

BET* 

1 . 3295E+32 

BET* 

1 .3321E+02 

BET* 

1 . 3295E+02 

BET* 

1.321 6  E + 3 2 

BET® 

1 . 3085E+O2 

BET* 

1 . 2903E+32 

BET® 

1 . 2669E+32 

BET* 

1 . 2336E+02 

BET* 

1 . 2054E  +  O2 

BET* 

1. 1 67  4E+02 

BET* 

1 .  1 248E+02 

BET* 

1 . 8777E+02 

BET* 

1 . 0264E+32 

BET* 

9. 7 1 03E  +  3 1 

BET* 

9. 1 191E+01 

BET* 

3 .  4914E+31 

BET* 

7. 8301 £+01 

BET* 

7. 1379E+01 

BET* 

6. 4176E+01 

BET* 

5 . 6723E+0 1 

BET* 

4 . 9O39E+0 1 

BET* 

4.  1  165E+3S 

BET* 

3. 3129E+01 

BET* 

2. 4962E+31 

BET* 

1  .  6  6  9  b  E  +  0 1 

BET* 

8. 3645E+00 

BET® 

1  .  3  3  2 1 E  -  0  7 

Pl*-4. 7455E-02 
PI*  1 . 5  0  4  0  E  -  y  2 
Pt*  S.133SE-02 
Ft*  1.511 7E-01 
PI*  2 . 2427E-01 
PI*  3 . 0034E-0 1 
PI®  3 . 7909E-0 1 
PI*  4.6020E-01 
Pi*  5 . 4 33b E -01 
PI*  6.2 S23E-01 
PI*  ?. 144SE-01 
PI*  8.01 7SE-0 1 
PI*  S.3977E-01 
PI*  9.7S10E-O1 
PI*  1.0664E+00 
PI®  1 . 1 544E+00 
PI*  1 . 241 7E+00 
Pl=  1 . 3230E+00 
PI*  1.41 2SE+00 
PI*  1.4960E+00 
PI®  1.577  IE +  00 
PI®  1 . b359E+0G 
P  1  =  1.7  319E+00 
PI®  1.8050E+0O 
PI*  1 . 3  7  4  9  E  +  0  0 
PI*  1.941 2E+00 
Pl=  2. O037E+00 
PI*  2 . 0621 £+00 
PI*  2.1 153E+80 
PI*  2.  l€b0£+08 
PI®  2 . 2109E+08 
PI®  2.251 1E+00 
PI®  2 . 23b 2 £+00 
PI*  2. 21EIE+00 
Pt®  2, 3408E+00 
PI®  2 •  3 b 0 0 £+00 
PI*  2. 3739E+80 
PI*  2 . 3322E+00 
PI*  2 . 3350E+00 


P2*  1.1731E-02 
P2*-l . 555 1 E- 02 
P2*-3.2355E-02 
F2*-3. 6541E-02 
P  2  *  -  2 .  b  0  7  0  E  ~  0  2 
P2=  9.6031E-04 
P2=  4. *2736-02 
P2®  1.11 3SE-0 1 
P2=  1 • 9754E-01 
P2*  3.O57bE-01 
P2=  4 . 3672E-0 1 
P2®  5.9030E-01 
P2*  7 • b303E-0 1 
P 2*  9.631OE-01 
P2®  1  . 1 904E+0G 
P2=  1 • 433SE+00 
P2=  1 . 897 1 E  +  00 
P2®  1 . 9785E+00 
P2®  2 • 2760E+00 
P2=  2 . 587 4 E  +  00 
P2*  2.9101E+8 0 
P2®  3.241  3E  +  8C 
P 2  =  3 . 5780E+  0 0 
P 2®  3 . 9 1 78E+0O 
P2®  4 . 255 1  E+00 
P2*  4. 53S3E+80 
P  2  *  4 . 9 1 4  9  E  +  0  0 
P2*  5. 2299E+00 
F2=  5 . 53336+00 
P2*  5.81 35E+00 
F2*  6.8759E+00 
c,2*  *.3147E*00 
P  2  ®  6  •  5  2  7  4  E  +  0  0 
,  7  1  1  8  E  ^  0  0 
P  2  *  8  .  $  £  5  2  E  +  0  0 
P  2  *  6 . 9  8  6  5  £  +  &  0 
P2®  7.0742E+00 
P2*  7 . 1 272E  +  33 
P2=  7 . 1 449E+0O 


r 


WnT »  1.1933E+O0 

DET* 

6.3016E+31 

PI*- 3. 1388E-81 

P2*  7.3443E-01 

WnT*  1.2566E+O0 

DET* 

6 . 2954E  +  G  1 

P 1 =-3 »  2397E-0 1 

F‘2*  3.  1  124E-01 

WnT  *  1 . 3 1 93E+00 

DET* 

6.2655E+01 

P 1 =  -3 , 3470E-0 1 

P'2*  3 . 9 1 34E-0  1 

WnT*  1  . 3323E+O0 

DET* 

6 . 2 1 34E+0 1 

P 1 =“3 . 4595E-0 1 

P2=  9 . 7402E-O1 

WnT  =  1 • 445 i  E+O0 

DET* 

6. 14O7E+01 

P 1 = - 3 , 5760E-01 

P"2*  1 . 0586E ^00 

WnT*  1 . 5080E+00 

DET* 

6 . 049QE+0 1 

PI— 3. 6956E-01 

P2*  1 . 1 443E+00 

WnT*  1 . 5703E+O0 

DET* 

5. 9397E+0 1 

F‘  1  =-3 . 8 i 73E-0 1 

P2=  1.2305E+00 

WnT*  1.6336E+00 

DET* 

5.8144E+01 

P  l  =  -  3 . 9  4  0  2  E  -  O  1 

P2*  1.31 66E+08 

WnT*  1.6965E+00 

DET* 

5 . 6744E+0 1 

P  l  =  -  4  . 0  6  3  6  E  -  0 1 

P2=  1.401 SE+00 

WnT*  1 . 7593E+O0 

DET* 

5, 5213E+01 

Pi— 4. 1867E-01 

P2=  1.48S7E+00 

WnT*  1.8221E+O0 

DET* 

5, 356  3 E  +  0 1 

P 1  * - 4 . 3090E-0 1 

P2*  1.S676E+00 

WnT*  1 . 8350E+O0 

DET* 

5. I810E+01 

P 1 =~4 . 4  2  9  9 E - 0 1 

P2*  1 . 64? 1E+00 

WnT*  1.9478E+80 

DET* 

4. 9964E+01 

Pi— 4. 5439E-01 

P2*  1 . 7236E+00 

WnT*  2.0106E+00 

DET* 

4. 3041E+01 

Pl=-4. 6657E-01 

P2*  1 . 7969E+O0 

WnT  *  2 . 0735E  +  O0 

DET* 

4 . 6050E+0 1 

PI— 4. 7799E-01 

P2=  1.3666E+0G 

WnT*  2 . 1 383E+00 

DET* 

4. 4005E+0 1 

Pi  — 4.3912E-01 

P 2  *  1 . 9  3  23  E  +  0 y 

WnT*  2.1991E+00 

DET* 

4.  191 5E  +-0 1 

P 1 * - 4 . 9993  E - 0 1 

p  2*  1 . 9938E  +  00 

WnT*  2.2619E+00 

DET* 

3. 9793E+01 

P 1 = - 5 . 1 04 1 E-0 1 

P2=  2.051 0E+G0 

WnT*  2 • 3243E+00 

DET* 

3 . 764-3E  +  0 1 

Pi— 5. 2054E-01 

P 2 *  2. 1036E  +  00 

W n T  *  2  m 3376E+00 

DET* 

* . 5439E+0 1 

P 1 = - 5 . 303  IE-01 

P2 *  2.151 &E  +  00 

WnT*  2. 4504E+W0 

DET* 

3. 3326E+01 

P  1  =  -  5 . 2  9  7  2 E  -  0 1 

P  2  =  2 .  1  9  4  9  E  ♦  0  0 

WnT*  2 . 5 1 3 3 E  +  0 0 

DET* 

3.  1  16SE+01 

P  1  —  5 . 4  8  7  6  E  -  0 1 

P2  =  2. 2 33 4 E  + 00 

WnT*  2 •  5  7  6  1  E  +  0  0 

DET* 

2 . 902  1 E  +  0 1 

Pi— 5. 5742E-01 

P2=  2 . 2 6 7 4 E  > v-uj 

WnT*  2.6339E+O0 

DET* 

2. 6S95E+01 

P  1  —  5  . 657OE-01 

P  2=  2.2966E  +  00 

WnT*  2. 7018E+00 

DET* 

2. 4795E+0 1 

P 1 = -5 . 7362E-0 1 

P2*  2 . 3 2 1 4 E  +  O 0 

WnT*  2. 7646E+GG 

DET* 

2.27 27E+01 

Pi— 5. 81 13E-81 

P2=  2 . 34 1 7E^00 

WnT*  2 . 3274E+O0 

DET* 

2 . 0699E+0 1 

P 1 = - 5 . 38  37E -0 1 

F'2*  2. 3578E  +  00 

WnT*  2.3903E+OO 

DET* 

1 . 87  X SE+01 

PI— 5.9521E-01 

F'2*  2 . 3 6 93 E  +  00 

WnT*  2 . 9 5 3 1 E  +  0 O 

DET* 

1 . 6731 E  + 0 1 

Pla-S.0i7lE-0i 

P2=  2. 3778E+0© 

WnT*  3.0159E+00 

DET* 

1 . 4  9 0  0  E  +  0  1 

P  1  *  -  6 . 0  7  3  7  E  -  0  1 

F‘2*  2.3822E  +  00 

WnT*  3. 0788E+8O 

DET* 

1 . 3077E+0 1 

P  1  =  -  6  .  1  3  7 1 E  -  0 1 

F’2*  2.3830E  +  00 

WnT*  3 . 1 4 16E+00 

DET* 

I  .  1  3 1 4  E  *♦*  0 1 

P 1 = - 6 . 1 923E-0 1 

F'2  *  2 . 3305 E +  O0 

COEFFICIENTS  OF 

THE 

characteristic  equation  with 

Pl»-2.37?QE-0l 

P2  =  5.307 

5E-01 

c  •  0  •  ♦  C  <  1  >  *  s  + 

Cc.  2 

)  «r  $  ■  •  '  2 

+  C  3 *  3  "v  3  + 

C<  0>»-4.?'?3301E-*'00 
C  <  l>»  I.316675E+01 
C<  2>*-l  .S40?99E+0l 
C'.  3  8.  1  18S82E+O0 

ROOTS: 

PE*L 

I 

mhG I nhPY 

D  P  M  F  i  n  0  P  R  7  I  0 

WnT 

1.61 1799E+00 

-1  . 

2J1697E-12 

-1.0 0 0 0 0 0 E  + 0 0 

4. 773508E-01 

3 , 3  9  3  3  S2E-0 1 

-4  . 

637730E-0 1 

4 . 9  S  4  6  4  5  E  -  0  1 

1 . 0  0  6  3  7  2  E  +  0  0 

3. 893332E-01 

4, 

6  37730E-01 

4.  ? 3 4 6 4 5 E - 0  1 

1  .  O  O  6  3720  O  O 

-  -  " 

4 

COEFFICIENTS  OF 

THE 

CHARACTERISTIC  EQUATION  WITH 

Pl--3.il  1 3E-01 

P2*  2. 9 6 0 0 E - 0 1 

C  1  0  ■  *  C  ■  1  *3  > 

C 's  2 

-  * 1  2 

♦  C  *  3  '  ■*  3  ‘  3  + 

C  G  1  s-4. 412458E+G0 
C  1  •  *  1  . 5376  l  9E  +  0  1 
r ■  2 «*-i . 940799E+01 

C-  3.113832E+00 


POOTS: 
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REAL 


8. 232520E-Q1 
7. 439712E-01 
3. 23252OE-01 


IMAGINARY 


-2. 297163E-8 1 
3. 331603E-12 
2. 297163E-01 


DAMPING  RATIO 


4. 99  758  3 E - 0 1 
1 . 000000E+00 
4. 997533E-01 


WnT 


3. 1 41594E-31 
2 . 9  5  7  5  3  0  E  -  0 1 
3. 141594E-31 


Y  o  u  h  au  e  selected  to  print  d  at  a  f  r  o  m  t-  h  e  p  ar  am  e  t  e  r  pi  an  e  p  r  o  g  r  am . 

This  program  placed  on  the  tape  the  values  of  P2  as  a  function  of 

PI.  You  will  be  given  the  maximum  and  minimum  values  of  these  variables 

but  you  can  select  any  value  you  wish  for  plotting  purposes. 

PI  is  e  qu  i  v  al  e  n  t  t  o  X  an  d  P  2  is  e  qu  i 1 1  a  1  e  n  t  t  o  Y . 


CASE  # 
1 


M  l  N  X 
-4. 29E-01 
-6. 19E-01 


MAX  X 
2. 3  3  E  +  0  0 
2.73E-01 


M  I N  V 
-3.S5E-02 
2. 74E-01 


MAX  V 
7. 14E+O0 
2. 38E+00 


f 


n--k< 


FIGURE  3  Parameter  Plane  Stability  Boundary  Sketch 


NOTE:  Single  cross-hatched  lines  are  z  *  +1  and  z  =  -1  boundaries 

Double  cross-hatched  curve  is  5  =  0  boundary. 

Encircled  number  represent  number  of  stable  roots 
in  indicated  region. 
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is  increased  and  shown  in  Fig.  5  to  facilitate  selection  of  a  robust 
design  point.  That  point  (underlined  in  the  printout)  is 


P 1  •  Kq  *  -0.32113 

and 

P2  *  K  ■  0.29601 . 

D .  Performance  of  System 

The  performance  of  the  candidate  design  is  evaluated  by 
applying  the  Cross-Multiplication  Technique  (Appendix  A)  to  the 
closed-loop  transfer  function  of  Eq.  ( 1 1 1  -  9 ) .  The  coefficients 
of  Eq.  ( 1 1 1  - 9 )  are  re-named  and  re-cast  as 

2 

.  £  n.  zJ 

ILL) —  =  - - -  .  ( 1 1 1  - 1  2 ) 

0c{z)  Z  h.  z^ 

j  =  0  J 


where 

h.  «  d.K  +  f  K  +  g  ( 1 1 1- 1 3 ) 

J  j  u  3  1  3 

and  the  numerator  coefficients  are  defined  in  Table  a. 


TABLE  A-  NUMERATOR  COEFFICIENTS 


0  Kgb  sin  uiT 

1  K^b  sin  o»T 

K2b  sin  ui T 


2 


3 

If  one  applies  Seltzer's  Cross-Multipl ication  Method  to  Eq.  (III-12), 
one  obtains  the  difference  equation, 

e(k+i)  =  k3[vc(|0  +  rijMk-1)  +  n0Mk-2) 

h  e(k-l)  "  h!^k-2)  -  %®(k-3)]  *  ( 1 1 1  -1 4) 

•  • 

where  9(k+l)  indicates  the  value  of  9(t)  at  the  instant  of  time, 
t  *  (k+l)T;  a  value  that  is  held  until  the  next  sampling  instant, 
t  =  ( k+2 )T;  etc . 

One  obvious  advantage  of  SAM  is  that  the  form  of  the  reference 
input,  © c ( t ) ,  need  not  be  specified  except  at  the  sampling  instants. 
Tables  5  and  6  indicate  the  missile  system's  response  to  unit 
step  and  ramp  inputs,  respectively.  They  are  shown  graphically  in 
Figures  6  and  7  . 


TABLE  5.  COMPUTER  PRINT-OUT  FOR  UNIT  STEP  RESPONSE 
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ORDER  OF  NUMERATOR  IS  2  ORDER  OF  DENOMINATOR  IS  3 

SAMPLE  INTERVAL  =  .1 


CALCULATIONS  MADE  ACCORDING  TO  FOLLOWING  EQUATION  DEFINITION 
C  <  Z  > •••■' R <  Z  )  =  NUMERATOR.  DENOMINATOR 

NUMERATOR  =  A  <  0  >  *Z'S0  +  A  <.  1  :•  *ZA  1  +  A<2J>*Z"2  .  .  .  A  <  N  >  * 
DENOMINATOR  «  Bm3>*Z'0  +  B<  1  ;*ZA1  +  B'2>*2'2  .  .  .  B<D>* 


NUMERATOR  COEFFICIENTS  A < 0 >  ,  A (1  )  ,  A < 2 )  ,  .  .  .  A ■; N > 
3. 7064  -3.4164  11542 


DENOMINATOR  COEFFICIENTS  B<0>  ,  B<  1  >  ,  B < 2 :• ,  .  .  .  B«  N> 
-4.41246  15.8762  -19.403  8.11839 


NUMBER  OF  STEPS  FOR  RESPONSE  =  38 
THE  INPUT  IS  A  UNIT  STEP  FUNCTION 


K 

TIME 

INPUT 

RESPONSE 

8 

0. 000E+00 

1 . 000E+00 

0. 000E+00 

1 

1 . 000E-01 

1 . 000E  +  0O 

-1 . 422E-82 

2 

2. 000E-01 

1 . 000E  +  00 

-4. 690E-01 

3 

3. 000E-01 

1 . 000E+00 

-1 . 072E+00 

4 

4. 00OE-01 

1 . 000E+00 

-1 . 631E+00 

5 

5. 00OE-01 

1.000E+0 0 

-2. 037E+00 

£ 

6. 0O0E-0 1 

1 . 0 00E+00 

-2. 241E+00 

7 

7. 000E-01 

1 . 080E  +  Q0 

-2. 238E+0O 

3 

S.000E-01 

1 . 000E+00 

-2. 053E+09 

q 

9. 000E-01 

1 . O00E+00 

-1 . 729E+O0 

10 

1 . 000E+00 

1 . 000E+00 

-1 . 313E+06 

1  1 

1  .  100E  +  00 

1 . 000E+00 

-8. 512E-01 

12 

1 . 200E+00 

1 . 080 E +08 

-3.861E-01 

1  3 

1 . 300E  +  0O 

1 . 000E+00 

4. 967E-02 

1  4 

1 . 4O0E+00 

1 . 0  0  0  E  +  0  0 

4 . 3 2  6  E  -  0 1 

15 

1 . 50OE  +  00 

1 . 0  0  0  E  +  O  0 

7. 486E-01 

16 

1 . 60OE  +  C0 

1  .  O00E  +  00 

9. 921 E -01 

17 

1  • 70OE  +  00 

1 . 000E  +  O0 

1  .  1 6  4  E + 0  0 

13 

1 . 3O0E  +  O0 

1 . 000E+00 

1 . 272E+O0 

19 

1 . 900E+00 

1  .  O08E  +  00 

1.324E+0O 

20 

2 . 000E+O0 

1 . 000E+00 

1 . 332 E +O0 

2  1 

2. 100E+00 

1 . 000E  +  00 

1 . 308E+0O 

22 

2 , 200E+00 

1 . 000E+00 

1 . 264E+0U 

2  3 

2 , 30OE+O0 

1 . 000E+00 

1 . 208E+00 

24 

2 , 400E+00 

1  .  O00E  +  00 

1 . 14SE+00 

25 

2. 500E+00 

1 . 000E  +  0O 

1 . 092E+00 

26 

2. 600E+00 

1 . 000E+80 

1 . 043E+00 

*■* 

4.  1 

2 . 700E+80 

1 . 000E+O0 

1 .0O3E  +  00 

23 

2. 900E+00 

1 .00OE  +  O0 

9. 732E-01 

29 

2. 900E+30 

1 . 000E+00 

9. 536E-01 

30 

3. 000E+00 

1 . 000E  +  00 

9. 429E-01 

IUMBER 

OF  STEPS  FOR 

RESPONSE  =  50 

’HE  INPUT  IS  A  UNIT 

STEP  FUNCTION 

K 

0 


TIME 

0. 000E+00 


INPUT 
1 . 0O0E+O0 


RESPONSE 
0. O00E+80 


r  j  r-j 


1 . 00OE-01 
2.O00E-01 
3. 000E-01 
4. 000E-01 
5. 000E-01 
6. 00OE-01 

r.  000E-01 
3-  000E-0 1 
9. 000E-01 
1 . 000E+00 
1 . I 00E+00 
1 . 200E+00 
1  .  -30OE+00 
l . 40OE+0 0 

1 . 7 00 E^ 00 
1 . 600E+G0 
1 . 700E+00 
1 . 300E+00 
1 . 9 0 0 E  4-  0  0 
2 . 000£>00 
2. 10OE+00 
2 , 2O0E+O 0 
2 . 300E+O0 
.  400E+00 
2. 300.E+00 
2 , €00E  +00 
2 . 7OOE+00 
2. 300E+00 
2. 9O0E+00 
3. O00E+00 
3. I00E>00 
?. 20OE+G0 
3. 300E+00 
3 . 4  0  O  E  +  0  0 
3,50 0 E  +  0 0 
3  ,  6  0  0  E  >  0  0 
3 . 7  0  0  E  >  O  0 
3 , 30  0 E  +  0 0 
3 . 9Q0E+00 
4. 000E+00 
4 . 1 0OE+00 
4 . 200E+O0 
4 .  300E  +  0  0 
4. 400E+00 
4 .  5OOE>00 
4 . 6  0  0  E  >  0  0 
4 . 70CE+00 
4, 300E+00 
4 . 9  0  0  E  +  0  0 
5. 000E+00 


1 . 000E  +  00 
1 , 000E+00 
1 , 000E+00 
1 . 0  0 0  E  +  0  0 
1 . 000E+00 
1 . 00GE+0O 
1 , 000E+00 
1 . 000E+00 
1 . 000E+00 
1 . 000E+00 
1 . 000Ef 00 
1 , 0 0 0 £  +■  0  0 
1 . 0  0  0  E  >  0  0 
1 . 000 E >00 
I . O0OE+0O 
1 . &00E+0O 

1 . 000E>00 
1 . 0 0  0 E  ♦  0  0 
2 , 000E+O0 
1 . 000E+00 
1 . 000E+00 
1  -  000E  +  00 
1 . 000E  +  00 
1 . 000E+00 
1 . 000E+00 
1 . 0  0  0  E  >  0  0 
1 . 0  0  0  E  >  0  0 
1  • 0O0E  +  0O 
1 . OOOE+00 
1 , 000E+00 
1 . 000 E >00 
1 . 000E+00 
1  •  000E>00 
1 . 0  0  O  £  >  0  0 
1 . 000E+00 
I . 000E+00 
1 • 000E>00 
1 . 000E>00 
1 . 000E+00 
: . 000E+OO 
1 . 000E>00 
1 . 00OE>00 
1 . 0O0E+0O 
I • 000E+00 
i . 000E+00 
2 . 000E+00 
1 . 000E>00 
2 . O00E+O0 
1 . 000E+00 
1 . 0  0  0  E  >  O  0 


-1.422E-02 

-4, 690E-01 
**1  *  07^E>0y 
**  1 . 631  £ >00 
-2. 037E+00 
■ 24 1  &  >00 
•  iw’S  E  >  0  0 
-2 . 053E+00 
“1 . 7  2 9  E  >  0  0 
*“1.31  3E  >00 
*3. 5 1 2  E - 0 1 
-3. 361 £-01 
4. 967E-0-2 
4. 326E-01 
7 . 4  3  6  E  -  G  l 
9. 9 21E-01 
1  .  1  6  4  E  >  0  0 
1 . 272E+0O 
1 . 3  2  4  E  >  O  o 
1 . 3  3  2  E  >  G  0 
1 . 303E+0O 
1  •  2  6  4  E  >  0  0 
1 . 2  0  3  E  >  0  O 
1  .  1  4  8  E  >  u  0 
1 . 092E>00 
1 . 043E>00 
1 . 0  0  3  E  >  O  0 
3  •  7  3  2  E  -  y  1 
'3. 536E-01 
9. 429E-01 
9. 39SE-01 
9, 425E-01 
9 , 4  9  3  E  ~  0 1 
9 .  5  8  4  E  -*  G  1 
9.685E-01 
9. 7  3  4  E - 0 1 
9.  3 7 5 E - 0 1 
9. 9  5 1 E  ~  0  1 
1  .  O0  1 E  > 0 0 
1 . 0  0  5  E  >  0  0 
1 .  O  0  3  E  >  0  0 
1 . 0  O  9  E  >  0  0 
l , O09E>00 
1 . 0  0  8  E  >  0  0 
1 . 007E>O0 
1 . 006E+0O 
1 . 004E>00 
1 . 0 0  2  E  0  0 
2 . 001E>00 
1 . 0 0 0 £>00 
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TABLE  6.  COMPUTER  PRINT-OUT  FOR  RAMP  RESPONSE 


ORDER  OF  NUMERATOR  IS  2  ORDER  OF  DENOMINATOR  IS 

SAMPLE  INTERVAL  *  .  1 


CALCULATIONS  MADE  ACCORDING  TO  FOLLOWING  EQUATION  DEFINITION 
C<Z>/R <Z)  *  NUMERATOR- DENOMINATOR 

NUMERATOR  =  fU0>«ZA0  +  A< 1>*2~1  +  A<2->*Z'-2  .  .  .  A<N>*Z-'N 
DENOMINATOR  -  B <0>*Z'0  +  B<i:>*Zy'l  +  B<2>*2-'2  .  .  .  B<D>*Z'D 


NUMERATOR  COEFFICIENTS  A  O  > ,  A  >;  1  1 ,  A  <■  2  >  ,  .  .  .  A  <  H  » 
3. 7064  -3.4164  -. 11542 


DENOMINATOR  COEFFICIENTS  B<  0  ■  .  B  •  1  '  ,  B  •  it ,  .  .  .  E  '  N  > 
-4.41246  15.3762  -19.408  8.11889 


NUMBER  OF  STEPS  FOR  RESPONSE  =  50 
INPUT  IS  RAMP 
SLOPE  IS  =  . 1 


K 

TIME 

INPUT 

RESPONSE 

e 

0.000E+00 

0. 000E+00 

0 . 0  0  0  E  +  0  0 

i 

1 . 00UE-01 

1 , 000E-02 

0 . 0  0  O  E  +  0  0 

2 . 000E-0 1 

2 . 000E-02 

-1 . 422E-04 

3 

3.000E-O1 

3. 000E-O2 

-4. 332E-03 

4 

4. 000E-01 

4 , 0  0  0  E  -  0  2 

-1 . 535E-02 

cr 

5. 00OE-01 

5 . 000E-02 

-3. 186E-02 

6. 000E'01 

6. 000E-O2 

-5. 223E-02 

7 

7. 000E-01 

7.000E-02 

-7. 4  £  4  E - 0  2 

Q 

3 . 000E-0 1 

8. 0O0E-02 

-9. 702E-02 

9 

9. 0  0  0  E - O 1 

9. 000E-82 

-1 . 175E-01 

10 

1 . 000E+00 

1 . 000E-01 

“1 .  3  4  3  E - 0 1 

1  1 

1 . 100E+00 

1  . 100E-01 

-1 . 480E-01 

12 

1 .200E+00 

1 . 200E-01 

- 1 . 5E5E-0 1 

13 

1 . 300E+00 

1 . 300E-01 

-1 . b  0  3  E - 0 1 

14 

1 . 40QE+00 

1 . 400E-01 

-1 . 598E-01 

15 

1 . 500E+00 

1 . 500E-01 

-1 . 555E-01 

1  6 

1 . 600E+00 

1 . 600E-01 

-1 . 48OE-01 

17 

1 . 700E+00 

1 . 700E-01 

-1 . 381E-01 

13 

1 . 300E+00 

1 . 800E-O1 

-1 . 2S5E-01 

19 

1 . 900E+00 

1 . 900E-0 1 

-1 . 137E-01 

20 

2. 000E+00 

2 . 0  0  8  E  -  0 1 

-1  . O05E-01 

21 

2. 100E+0O 

2. 100E-O1 

-8. 7 1 9E-02 

•i  ■*, 

i.  4. 

2 . 200E+0O 

2. 200E-01 

-7.4UE-02 

•t  *. 

4. 

2 . 300E+00 

2. 300E-01 

-6  •  1 4  7  E - 0  2 

24 

2. 400E+00 

2 . 400E-O 1 

-4. 940E-02 

25 

2. 5 0 0 E 0 0 

2.500E-W1 

-3. 791E-02 

29 


2. 60QE+00 

2  •  60UE-kJ  1 

-2. 

27 

2 . 700E+00 

2. 700E-01 

-1 . 656E-02 

23 

2 . 800E+00 

2. S00E-01 

-6. 530E-03 

29 

2 . 900E+00 

2. 900E-01 

3. 202E-03 

30 

3. 000E+00 

3. 000E-01 

1 . 274E-02 

31 

3. 100E+00 

3. 100E-01 

2.21 7E-02 

o  ■“> 

3. 20QE+00 

3. 200E-01 

3. 157E-02 

33 

3. 3  0  0  E  +  0  0 

3.300E-01 

4. 099E-02 

34 

3. 400E+00 

3. 400E-01 

5. 043E-02 

35 

3 • 300E+00 

3. 500E-01 

6. 007E-02 

36 

3 • 600E+00 

■3. 600E-01 

6. 975E-02 

37 

3. 700E+00 

3. 700E-01 

7. 954E-02 

33 

3 . 300E+00 

3. 80OE-01 

8. 941E-02 

39 

3. 980E+00 

3. 900E-01 

9. 936E-02 

40 

4 . 000E  +  00 

4. 000E-01 

1 • 094E-0 1 

41 

4. 100E+00 

4. 100E-01 

1 . 194E-01 

42 

4. 200E+00 

4. 200E-01 

1.295E-01 

43 

4 . 300E+00 

4. 300E-0 1 

1 . 396E-01 

44 

4. 400E+00 

4. 4O0E-01 

I . 497E-01 

45 

4 . 500E+00 

4. 500E-01 

1 . 598E-01 

46 

4. 60QE+00 

4. 6 O 0 E - 0 1 

1 . 698E-01 

47 

4 . 70UE+00 

4. 700E-01 

1 . 799E-0 1 

48 

4. 300E+00 

4. 800E-01 

1 . 8  9  9 E  -  0 1 

4  9 

4 . 900E+00 

4. 900E-01 

2. 0 0 0 £  —  0 1 

50 

5. 00OE+00 

5. 000E-01 

2. 100E-01 

INPUT 

SLOPE 


IS  RAMP 
IS  =  .  1 


K 

TltlE 

INPUT 

RESPONSE 

0 

0. 000E+00 

0  .  0  0  0  E  +  0  0 

0 . 000E+00 

1 

1 . 000E-0 1 

1 . 000E-02 

0 . 0  0  8  E  +  0  0 

2 

2 . 000E-0 1 

2 . 0  0  0  E  -  0  2 

-1 . 422E-04 

j! 

3. 000E-01 

3 . 00OE-02 

-4. 832E-83 

4 

4 . 000E-0 1 

4 . 000E-02 

-1  . 555E-02 

Cj 

5. 000E-01 

5 . 0  0  0  E  -  0  2 

-3. 136E-02 

6 

6. 000E-01 

6.000E-02 

-5. 223E-02 

7 

7. 00OE-01 

7. 000E-02 

-7. 4  6  4  E - O  2 

3 

8. 000E-01 

3. 000E-O2 

-9. 7O2E-02 

9. 000E-01 

9 . 0O0E-82 

-1 . 175E-01 

10 

1 . 000E+00 

1 . 0O0E-0 1 

-1 . 34SE-01 

1  1 

1 . 1 00E+00 

1  .  1  00E- 0 1 

-1 . 4  3  0 E - 0 1 

12 

1 . 200E+00 

1 . 200E-01 

-1 . 5 6  5  E - 0 1 

13 

1 . 300E+00 

1  .  3  0  0  E  *  0  1 

-1 . 6O3E-01 

1  4 

1 . 400E+00 

1 . 400E-O1 

-1 . 598E-01 

15 

1 .500E+00 

1 . 500E-01 

-1 . 555E-01 

16 

1 . 600E+00 

1 . 600E-OI 

-t . 430E-01 

17 

1 . 700E+00 

1 . 70OE-01 

-1 . 381E-01 

18 

1 . 800E+00 

1 .300E-01 

-1 . 265E-01 

19 

1 . 900E+00 

1 . 9  0  0  E  -  O  1 

-1 . 137E-01 

20 

2. 000E+00 

2. 000E-01 

- 1 . 005E-01 

21 

2. 100E+00 

2. 1O0E-01 

-8. 719E-02 

->  -i 

4.  A. 

2. 2 O 0 E *»■  0 O 

2  •  2  0  0  E  -  O  I 

-7. 41  IE-02 

23 

2. 300E+80 

2. 3O0E-0 1 

-6. 147E-02 

24 

2. 4O0E+00 

2. 40OE-01 

-4. 940E-02 

25 

2 . 500E+00 

2.500E-O1 

-3.791E-02 

26 

2 . 600E+O0 

2.60OE-01 

-2. 699E-02 

&  i' 

2. 700E+00 

2. 700E-01 

-1 . 656E-02 

23 

2 . 300E+00 

2. 300E-01 

- 6 . 530E-03 

23 

2. 900E+00 

2. 900E-01 

3 4 202E-03 

30 

3. 000E+00 

3. O00E-01 

1 . 274E-02 

31 

3 . 1 00E+00 

3. 100E-01 

2.217E-02 

32 

3 . 2O0E+00 

3.200E-01 

3. 157E-02 

3  3 

3. 300E+00 

3. 300E-01 

4. 899E-02 

34 

3. 400E+O0 

3. 400E-0 1 

5 . 043E-02 

35 

3. 50OE+00 

3 • 5O0E-01 

6 . 00  7E-02 

36 

3.600E+O0 

3 . 600E-0 1 

6. 975E-02 

3? 

3 . 700E+00 

3. 70OE-01 

7. 954E-02 

33 

3. 300E+O0 

3. S00E-01 

3. 941E-02 

39 

3. 90QE+00 

3. 900E-01 

9. 936E-02 

40 

4. 000E+00 

4 . 000E-01 

1 . 094E-01 

41 

4, 100E+0O 

4 . 1O0E-01 

1  .  194E-01 

42 

4. 200E+00 

4.20  0  E - 0 1 

1 . 295E-01 

4  3 

4 . 300E+0O 

4. 300E-01 

1 . 3  9  6  E  -  0  1 

44 

4. 40QE+00 

4. 400E-O1 

1  •  4  9  7  E  -  0 1 

45 

4 . 500E+00 

4. 500E-01 

1 . 593E-01 

46 

4. 600E+00 

4 . 600E 'O 1 

1 . 693E-01 

47 

4, 700E+00 

4 • 700E~0 1 

1 . 7  9  9  E  -  0 1 

43 

4  *  300E  +  00 

4. 800E-01 

1 . 3  9  9 E  -  0 1 

49 

4  *  900E  +  00 

4. 900E-01 

2. 000E-01 

50 

5 , 000E+00 

5.000E-01 

2. 10OE-01 

51 

5, 100E+00 

5.  100E-01 

2. 200E-01 

52 

5 , 200E+00 

5. 200E-01 

2. 300E-01 

C  -1 

5 . 3  0  0  E  +  O  0 

5. 300E-O1 

2. 400E-01 

54 

5. 400E+00 

5 . 400E-0 1 

2. 499E-01 

55 

5 . 500E+00 

5. 500E-01 

2. 599E-01 

56 

5. 600E+0U 

5. 600E-O1 

2. 699E-01 

57 

5. 700E+00 

5 . 700E-0 1 

2. 799E-01 

53 

5. 3O0E+00 

5. 300E-01 

2. 899E-01 

59 

5. 900E+00 

5. 900E-01 

2. 999E-01 

60 

6. 00QE+00 

6. 0O0E-01 

3.099E-01 

6  1 

6. 100E+00 

6. 100E-O1 

3. 19SE-01 

62 

6. 200E+00 

6. 200E-01 

3. 298E-01 

6  3 

6. 300E+00 

6. 3 0 0 E - O  1 

3. 398E-01 

6  4 

6. 400E+00 

t  6.400E-01 

3. 498E-01 

65 

6 . 500E+00 

6. 50OE-01 

3. 59SE-01 

66 

6 . 600E+00 

6. 6 0 0 E - 0  1 

3 .  $  9  3  E  -  O  1 

67 

6. 708E+00 

6.  700E-01 

3. 798E-01 

63 

6 . 300E+00 

6. 3O0E-01 

3. 393E-01 

69 

6. 900E+00 

6. 90OE *0 1 

3. 99SE-01 

70 

7.000E+00 

7. 000E-01 

4. 098E-01 

71 

7. 10OE+0O 

7. 100E-O1 

4. 19SE-01 

7  2 

7. 200E+00 

7. 200E-01 

4. 29SE-01 

73 

7. 300E+00 

7. 300E-01 

4. 39SE-0 1 

7  4 

7 , 4O0E+O0 

7. 400E-01 

4. 49SE-01 

—  e- 

-  -j 

7 . 500E+00 

7. 500E-01 

4. 59SE-01 

76 

7 . 600E+00 

7 . 6  ©  0  E  -  O  1 

4. 69SE-01 

7  7 

7. 70OE+0O 

7. 700E-01 

4. 79SE-01 

70 

7. 300E+00 

7. 300E-O1 

4.898E-01 

7  9 

7. 900E+00 

7 • 900E-0 1 

4. 993E-01 

30 

9. O00E+00 

3. 000E-01 

5. 093E-01 

31 

3. 1O0E+00 

3. 100E-01 

5. 198E-01 

8 . 200E+00 
3. 300E+00 
8. 400E+00 
8. 500E+08 
3. 600E+80 
8. 700E+00 
3 . 800E+00 
8. 308E+0O 
9. 000E+00 
9. 1 00E+00 
9 . 20OE+00 
9. 3O0E+00 
9. 40OE+00 
9 . 500E+00 
9 . 60OE+80 
9. 700E+00 
9 . S80E+00 
9. 980E+00 
1.000E+01 


8. "200E-0r 

8. 300E-01 
3. 400E-01 
8.500E-01 
3. 680E-01 
8 . 700E-0 1 
8. 800E-01 
3. 900E-01 
9. 000E-01 

9.  100E-01 
9. 200E-01 
9. 300E-O1 
9. 400E-01 
9. 500E-01 
9. 800E-O1 
9. 700E-01 
9. 800E-01 
9. 900E-01 
1 . 000E+00 


5T298E-0 1~ 
5. 398E-01 
5. 493E-01 
5. 598E-01 
5.69SE-01 
5. 798E-81 
5. 898E-01 
5. 99SE-01 
6. 093E-01 
6.  198E-81 
8. 293E-01 
8. 398E-01 
8. 498E-01 
8.593E-01 
8. 898E-01 
8 . 798E-0 1 
8.39SE-01 
6. 997E-01 
7.097E-01 
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SECTION  IV.  CONCLUSIONS 

The  planar  equations  of  motion  have  been  developed  for  the 
FAKMS  missile  system.  They  were  linearized  to  permit  stability 
and  preliminary  autopilot  design  to  proceed.  A  technique  for 
performing  this  preliminary  design  has  been  developed  and  applied. 
In  the  actual  case,  the  design  would  be  repeated  until  a  near- 
optimal  design  was  achieved.  It  would  then  be  repeated  for 
several  linearized  operating  points  along  the  path  of  a  nominal 
trajectory.  Then  the  system  would  be  "flown"  on  a  simulation 
with  numerical  values  for  autopilot  coefficients  changed  during 
flight  if  necessary.  Finally,  the  procedure  would  be  repeated 
for  three-dimensional  eauations  of  motion. 
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APPENDIX  A.  DETERMINATION  OF  DIGITAL  CONTROL  SYSTEM  RESPONSE 
BY  CROSS-MULTIPLICATION 

The  following  material  was  developed  under  US  Army  Missile 
Command  Contract  DAAK40-78-C-0226  and  completed  under  this 
contract  (DAAH01 -81-M-A1 1 8) .  It  was  accepted  and  presented 
at  the  IEEE  Region  3  " SOUTHEASTCON  ’81"  Conference  held  in 
Huntsville,  Alabama  on  5-8  April  1981. 


DETERMINATION  OF  DIGITAL  CONTROL  SYSTEM  RESPONSE  BY  CROSS-MULTIPLICATION 


SHERMAN  M.  SELTZER 


Control  Dynamics  Comnany 
Huntsville,  Alabama 


ABSTRACT 


II.  RESPONSE  AT  SAMPLING  INSTANTS 


This  paper  describes  the  Cross-Multiplica¬ 
tion  technique  for  obtaining  the  response  of  a 
digital  control  system  from  its  closed-loop 
transfer  function.  When  the  latter  is  expressed 
in  the  complex  z-domain,  the  response  may  be 
obtained  at  the  sampling  instants.  If  it  is 
desired  to  know  the  response  between  the  sampling 
instants,  the  Submultiple  Method  may  be  adapted 
to  the  Cross-Mul tipi ication  technique.  When 
the  transfer  function  is  expressed  in  the  Mod¬ 
ified  z-transfonm  domain,  the  response  is 
obtained  at  any  point  within  each  sampling 
instant.  The  three  aDDroaches  are  described 
and  applied  tc  an  example. 


r.  INTRODUCTION 

This  paper  describes  a  technique  for  obtaining 
the  response  of  a  digitally  controlled  system. 

It  is  assumed  that  the  closed-loop  transfer  fun¬ 
ction  is  available  in  the  z-or  modified  2-trans¬ 
form  domain.  The  numerator  and  denominator  of 
each  side  of  the  transfer  function  are  cross- 
multiolied.  The  Real  Translation  Therorem  is 
then  applied  tc  the  result,  yielding  a  difference 
equation  in  the  time-domain.  This  may  be  solved 
for  the  system  response  in  terms  of  the  refer¬ 
ence  (or  other)  innut(s)  to  the  system  as  well 
as  in  terms  of  system  state  initial  conditions. 

Two  different  modifications  to  the  basic  tech¬ 
nique  are  described:  one  usina  the  submultiple 
method  and  one  usinq  the  modified  z-transform 
technique.  These  are  applied  when  it  is  desired 
to  determine  intra-sampl inq  responses  of  the 
system.  All  three  techniques  are  applied  to  a 
single  example.  A  summary  of  the  techniques 
and  their  application  is  provided  at  the  conclu¬ 
sion  this  note. 


It  is  assumed  that  a  given  digital  or  sampled- 
data  system  can  be  described  by  a  closed- loop 
transfer  function  that  relates  the  controlled  out¬ 
put  of  the  system  to  the  reference  inout.  If 
there  is  more  than  one  input,  the  technique  can 
also  be  applied  to  the  resulting  sum  of  closed- 
loop  transfer  functions  relating  the  controlled 
output  to  each  of  the  inputs.  Although  this  note 
refers  only  to  a  single  controlled  output,  the 
technique  can  be  applied  to  find  any  system  state 
if  it  is  related  to  the  inputs  to  the  system  in 
the  z-domain. 


It  is  assumed  that  the  state  whose  response 
is  desired  is  denoted  in  the  z-domain  as  C(z)  and 
the  reference  input  as  R(z).  The  relationship 
between  C(z)  and  R(z)  usually  can  be  expressed  as 
a  closed-loop  transfer  function  (or  several  such 
transfer  functions)  which  is  a  ratio  of  two  poly- 


1  ,e. , 


4“0  J  _ 

N 


(U 


where  coefficients  a.  and  b^  represent  the  system 
parameters.  The  procedure  Tor  findinq  the  re¬ 
sponse  at  sampling  instants--by  the  cross-multi¬ 
plication  method--consi sts  of  three  steps.  First 
one  cross-mul tiDl ies  the  numerators  with  the  demon- 
inators  of  Equation  (1).  Then  each  term  of  the 
resulting  expression  is  divided  by  bNz".  This 
result  is  solved  for  the  C(z)  term  that  is  not 
multiplied  by  a  nonzero  power  of  z,  i.e. 


Finally,  the  Real  Translation  Theorem  is  applied 
to  Eauation  (2),  The  resulting  values  of  c(nT) 
at  each  sampling  instant,  nT,  become 
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The  advantages  of  the  form  of  Equation  (3)  are 
three- fold:  (1)  the  value  of  c(nT),  for  any  t=nT, 
may  be  obtained  for  any  form  of  r(t),  whether  or 
not  it  is  "z-transformable;M  (2)  the  expression 
for  c(nT)  does  not  have  to  be  recalculated  every 
time  the  form  of  r(t)  changes,  as  is  the  case  when 
the  response  is  determined  by  the  partial  fraction, 
power  series,  or  inversion  formula  methods;  and 
(3)  the  form  of  the  expression  for  c(nT)  permits 
the  inclusion  of  initial  conditions,  such  as 
c(0),  if  they  exist. 

III.  RESPONSE  BETWEEN  SAMPLING  INSTANTS: 

SUBMULTIPLE  METHOD 

If  it  is  desired  to  find  the  intra-sampling 
response  of  the  same  type  sampled-data  system 
described  in  Section  II,  it  may  be  accomplished 
by  applying  the  submultiple  method.'  Briefly, 
let  c(mT/n)  represent  the  value  of  the  response 
c ( t )  at  the  instant,  t=mT/r,  wbe^e  n-1  represents 
the  number  of  intrasampling  responses  desired 
(n  is  an  integer  with  value  greater  than  unity). 

If  m  is  an  integer,  the  sampling  period  within 
which  the  submultiples  are  to  be  determined  is 
denoted  as  mT.  The  z-transform  of  c(mT/'n)  may 
be  found  from  the  ordinary  z-transform  in  the 
following  manner.  In  essence,  c(mT/n)  is  the 
ouput  of  a  fictitious  sampler  which  samples  n 
times  as  fast  as  the  real  sampler.  The  z-trans¬ 
form  of  that  output  is  defined  as 


method.?  In  this  method  it  is  necessary  first 
to  determine  the  modified  z-transform  equivalent 
of  the  closed-loop  transfer  function  of  Equation 
(1),  l.e. 


where  a-  and  b.  may  be  functions  of  m.  The  cross- 
multi  pncationJtechnique  is  similar  to  those  two 
previously  described.  Again,  the  numerators  and 
denominators  of  Equation  (71  are  cross-  multiplied. 
Each  term  is  divided  by  bNzN,  and  the  resulting 
equation  is  solved  for  the  C{z,m)  term  that  Is  not 
multiplied  by  a  non-zero  power  of  z.  The  inverse 
z-transform  of  each  term  is  determined,  yielding 

c  [(n-l)T,mj  =  c  £<n- {&) 

1  C  N 

=  Cy  0-N+J>T] 

1  j*°  N-1 

’  r.b.  c  f  (N-l+m-N+k )T| 

k=0  L  J 

where  are  integers. 


^,^c(nT/m)|  d.  Cfz)^  *  C(z)| 
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where  z  =  zi/  and  T  =  T/n.  Now  the  closed-loop 
expression  of  Equation  (I)  may  be  altered  to  read 
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The  submultiple  modification  to  the  basic 
method  also  consists  of  three  steps  which  are 
analogous  te  the  steps  of  Section  II.  The  result 
is  the  expression 

c(mT/n)  -  1  fM  (6) 


“  TH.  bjC  [(m-N+j )T/i»j  1 
j=0  J 


where  all  values  of  rC*]  are  equal  to  zero  except 
tor  those  values  of  r[..]at  integral  multiples  of 
T,  at  which  time  r(pT)  eauals  the  value  of  r(t)  at 
t=pT;  j ,k ,M,m,N,n,p  are  integers. 

IV.  RESPONSE  BETWEEN  SAMPLING  INSTANTS: 

MODIFIED  Z-TRANSFORM  METHOD 


If  it  is  desired  to  find  the  intra-sampl inq 
response  of  the  same  type  of  sarr.pled-data  system 
described  in  Section  II,  it  also  may  be  accom¬ 
plished  by  applying  the  modified  z-transform 


V.  EXAMPLE 


Given  the  closed-loop  sampled-data  system 
shown  in  the  figure,  where  G( s)=l/s (s+I )  and 
H(s)sl.  The  closed-loop  transfer  function  is  fourf 
to  be 

C(z)  _  (1  -  e~T)  z _ 

R(z)  Z?  -  2e~TZ  +  e*T  (9) 


Applying  the  three  steps  prescribed  in  Section  II, 
one  may  obtain 

c(nT)  -  (l-e‘T)rC(n-l)T]  +  2e‘Tc  [(n-l)fj 


-  e"Tc  [(n-2)T^  . 


Now  values  for  the  response  may  be  found  at  the 
sampling  instants  by  lettirg  n  progress i vely  take 
on  values  of  0,1,2,  .  .  . 

If  it  is  desired  to  find  the  response  between 
sampling  instants  by  the  submultiple  method,  it 
may  be  shown  that 


C(z)n  =  G(z)pE(z)  *  G(z)nR{z)/  [l+HG(zj]  ,  (11) 


where  G{z)|jis  found  by  the  means  shown  in  Equation 
(4).  Application  of  the  procedures  of  Section  III 
yields  the  time-domain  difference  equation  for 
eft)  at  the  instant  of  time,  t=mT/n, where  m  is  any 
desired  integer.  A  fictitious  time  del ay9 e" ( 1 -m)T 
may  be  placed  following  the  forward  loop  gain, 
G(s),  and  a  fictitious  time  advance  (of  the  same 
magnitude)  may  be  place  just  before  the  feedback 
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loop  gain  H($).  The  closed-loop  transfer  function 
that  results  from  these  two  elements  being  added 
is 

(12) 

C(z,m)  =  G(z,ml  -  (l-e-mT)z  +  (e"mT-e-T)  . 

Ru)  1+  HGtz)  Z2  -  2e-T l  +  e-T 

It  may  be  observed  that  Equation  (12)  is  in  the 
same  form  as  Equation  (7).  Hence  the  prescribed 
operations  of  Section  IV  may  be  applied,  resulting 
in  the  difference  equation 

c  [(n-l+m)r]  =  (l-e"mT)  r[(n-l)T]  +  (e'mT-e“T)* 

•  r  [(n-2)T^J  +  2e~^  c  J(n-2+m)lJ 

•c  |(n-3+m)T]  .  (13) 

During  the  nth  sampling  period,  i.e.  (n-1)  T£ 
t-  nt,  where  n  is  an  integer,  one  lets  m 
assume  values  between  zero  and  one  to  find  the 
intrasampling  values  of  c(T). 

VI.  CONCLUSIONS 

Several  analytical  techniques  for  obtaining 
the  response  of  a  digital  control  system  have 
been  described.  They  are  based  on  a  sinqle 
principle:  cross-multiplication  followed 

by  applications  of  the  real  translation  theorems. 
Each  is  applied  to  a  single  example.  As  a 
starting  point  for  application  of  each  of  the 
techniques,  it  is  required  that  the  dynamics  of 
the  digital  control  system  be  described  in  the 
z-  or  modified  z-domain. 

The  advantages  of  the  three  techniques 
over  extant  classical  methods  are: 

The  response  may  be  obtained  for  any  deter¬ 
ministic  reference  input  into  the  system  as  long 
as  its  value  is  known  at  the  sampling  instants. 

It  need  not  be  described  by  a  differential 
equation,  and  the  z-transform  for  a  specific 
reference  input  need  not  be  determined  before 
obtaining,  an  expression  for  the  response. 

Initial  conditions  and  instantaneous  changes 
can  be  acconuiodated  readily  by  the  equations 
obtained  through  use  of  the  cross-multiplication 
methods . 

The  difference  equations  obtained  are  part¬ 
icularly  amenable  to  programming  on  a  desktop 
calculator  or  digital  computer, 

A  detailed  knowledge  of  the  theory  underlying 
digital  or  sampled-data  control  systems  is  not 
required  (although  it  certainly  is  helpful)  by 
the  analyst  in  order  to  apply  the  recipes 
described  herein. 
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Fiqure.  Closed-loop  sampled-data  system 
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APPENDIX  B.  A  VALIDATED  METHODOLOGY  FOR  ACCURATELY  PREDICTING 
MISSILE  FLIGHT  PERFORMANCE 

The  following  material  was  developed  in  part  and  completed 
under  this  contract  ( DAAH01 -81 -MAI  1 8) .  It  was  accepted  and 
presented  at  the  American  Astronautical  Society  Annual  Rocky 
Mountain  Guidance  and  Control  Conference  at  Keystone,  Colorado 
on  31  January  -  4  February  1981. 
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A  VALIDATED  METHODOLOGY  FOR  ACCURATELY  PREDICTING 
MISSILE  FLIGHT  PERFORMANCE 

Harold  L.  Past rick*  and  Sherman  M.  Seltzer t 

A  comprehensive  hardware- in -the-loop  (HWIL)  simulation 
has  been  developed,  operational  in  real-time,  to  per¬ 
form  guidance  and  control  system  design  verification 
and  autopilot  and  detector  optimization  for  a  laser 
semiactive  terminal  homing  weapon  system.  Emphasis 
was  placed  on  the  development  of  time-critical 
aspects  of  the  simulation  including  analog-digital 
computer  partitioning,  aerodynamic  moment  updating,  \ 
and  target-vehicle  geometry  modeling.  Validity  of  the 
simulation  was  demonstrated  by  comparing  results 
obtained  from  the  all-digital  6-DOF  and  hybrid  analog- 
digital  simulations.  The  simulation  was  proven  an 
invaluable  evaluation  and  optimization  tool  in  efforts 
to  enhance  weapon  system  performance,  since  simulation 
results  were  indicated  by  actual  flight  tests. 


INTRODUCTION 

It  is  a  widely  accepted  principle  in  the  missile  development  community 
that  the  ultimate  figure  of  merit  for  a  missile  is  its  performance  in 
flight.  Typically,  a  significant  amount  of  dollars  is  dedicated  to  the 
flight  test  phase  of  the  development  program  in  order  to  assess  its 
strengths  and  shortcomings  in  scenarios  ranging  from  simplistic  (advanced 
development  concepts)  to  battlefield  conditions  [Operational  Testing  (OT) 
under  .n  emulated  realistic  environment] .  To  the  extent  that  key  missile 
parameters  are  monitored,  the  principle  of  flight  testing  remains  as  the 
soundest  indicator  of  missile  system  performance. 

However,  it  is  not  unusual  to  have  a  missile  fail  in  a  flight  test.  Thus 
it  would  yield  little  or  no  data  for  reasons  such  as  telemetry  failure, 
camera  coverage  inadequacies  or  planned  nonmonitoring  of  parameters  as  is 
the  condition  for  OT  flights.  A  variety  of  reasons  such  as  these  and 
others,  to  yield  flight  equivalent  data,  motivates  the  search  for  an 
alternative  to  flight  testing.  At  the  US  Army  Missile  Command  (MICOM) , 
that  alternative  is  widely  accepted  under  the  name  of  hardware-in-the- 
loop  (HWIL)  simulation  testing.^  In  this  mode  of  simulation,  actual 


♦Chief,  Guidance  and  Control  Analysis,  Guidance  and  Control  Directorate, 
US  Array  Missile  Laboratory,  US  Army  Missile  Command,  Redstone  Arsenal, 
Alabama  35898. 

^Consulting  Engineer,  Control  Dynamics  Company,  221  East  Side  Square, 
Suite  IB,  Huntsville,  Alabama  35801.  •  -  ~ 


(1) 


41 


missile  flight  hardware  such  as  autopilot’,  actuator "seeker /inertial 
components  and  flight  computer  are  included  in  a  simulation  that  is  closed 
around  a  digital /analog  computer  that  has  been  coded  with  appropriate 
models  to  yield  an  emulated  flight  in  real-time. ^ 

The  intent  of  this  paper  is  to  summarize  the  concept  of  HWIL  simulation 
testing  as  performed  within- several  elements’  of  MICOM  amTliisplay  the 
correlation  of  those  tests  (which  are  nominally  nondestructive)  to  flight 
tests.  In  particular,  data  will  be  provided  to  show  that  hundreds  of 
preflight  runs  in  a  well-designed  HWIL  test  not  only  predict  flight  per¬ 
formance  to  within  3  sigma  tolerances,  but  actual  HWIL  simulation  tested 
hardware  performs  better  in  flight  because  of  the  advantage  of  having 
been  ”tweaked-to-per feet ion”  in  the  simulation  laboratory. ^  This  does 
not  imply  that  perfection  has  been  achieved  via  simulation.  To  the  con¬ 
trary,  as  the  example  of  Mlessons  learned”  has  shown. *  Post-flight 
analysis  of  two  early  guided  projectiles  indicated  aerodynamic  maneuvera¬ 
bility  was  less  than  that  predicted  by  wind  tunnel  tests  and,  of  course, 
those  aero  data  were  incorporated  in  the  simulation  models.  The  cause, 
effect,  and  solution  will  be  discussed.  Nevertheless,  savings  in  the 
millions  of  dollars  have  been  demonstrated  by  virtue  of  using  the  HWIL 
simulation  testing  as  a  precursor  to  flight  test. 

SIMULATION  EVOLUTION 

The  Guidance  and  Control  Directorate  of  the  US  Army  Missile  Laboratory, 
MICOM,  has  been  engaged  for  a  number  of  years  In  the  area  of  semiactive 
laser  terminal  homing  guidance.  In  essence,  .this_type  of  guidance  con¬ 
sists  of  designating  the  target  by  a  ground-based  laser  beam  and  using  a 
missile-mounted  laser  receiver/detector,  in  the  terminal  phase  of  the 
missile’s  flight,  to  provide  guidance  and  control  signals  by  tracking  the 
laser  energy  reflected  by  the  target. 

As  part  of  the  advanced  development  effort  applied  to  this  research  area, 
an  extensive  program  of  simulation  studies  was  undertaken  to  identify 
potential  problem  areas  and  to  indicate  approaches  to  their  solution. 

The  initial  approach  consisted  of  a  full  six  degree-of-f reedom  (6  DOF) 
digital  computer-based  simulation  which  focused  on  the  effects  upon 
missile  performance  of  variations  in  autopilot  design  during  the  evolu¬ 
tion  of  the  missile’s  configuration. 

The  next  stage  in  the  simulation  activity  took  place  over  a  longer  period 
and  was  related  to  the  choice  of  the  guidance  and  control  system  and 
optimization  of  the  guidance  law  and  autopilot  mechanization.  Amendments 
and  improvements  were  made  to  the  simulation  as  new  data  (aerodynamic 
coefficients,  hardware  parameters,  etc.)  became  available.  It  was, 
however,  realized  that  the  models  employed  were  idealized  and  suffered 
limitations  in  their  representation  of  the  real  world.  In  an  attempt  to 
overcome  some  of  the  more  obvious  limitations,  test  results  from  an  early 
version  of  the  laser  seeker  were  included  in  the  simulation  in  place  of 
the  idealized  seeker  model;  the  test  results  involved  the  introduction 
of  Gaussian  noise,  thus  requiring  a  Monte  Carlo  approach  to  obtaining 
results . 
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Further  efforts  to  investigate  and  validate  the  guidance  and  control 
system  were  directed  to  a  time  frame  immediately  after  the  design  stage 
was  complete,  when  prototype  hardware  had  been  produced  but  prior  to 
initial  flight  tests  of  the  missile.  The  method  selected  was  a  natural 
extension  of  the  earlier  simulation  work  and  consisted  of  a  hybrid  digital- 
analog  real-time  computer^  simulation  empl oying_t he  .guidance  and  control 
HWIL.  This  was  considered  an  important  part  of  the  development  phase 
because,  although  the  hardware  had  been  manufactured,  it  was  still  possi¬ 
ble  to  vary  a  number  of  the  design  parameters,  and  in  addition  optical- 
glmbal  coupling  in  the  seeker  was  recognized  as  a  problem  which  was 
difficult  to  handle  theoretically.  Furthermore,  closed  loop  operation  in 
a  realistic  simulation  offered  the  opportunity  of  studying  the  dynamic 
stability  and  performance  of  the  hardware  system  first  hand,  with  the 
hope  of  reducing  the  length  and  cost  of  the  flight  testing  cycle. 

HARDWARE- IN-THE -LOOP  SIMULATION  OVERVIEW 


The  hybrid  simulation  facility  consists  of  the  major  components  shown  in 
Fig.  1.  The  interfaces  contain  the  hybrid  input /output  and  communica¬ 
tion  devices  which  include  analog-digital  converters  (ADC’s),  multiplying 
digital-analog  converters  (MDAC’s),  discrete  digital  input /output  lines, 
interrupt  lines,  and  real-time  clocks. 


Fig.  1  Block  Diagram  of  hardware-in-the-loop  (HWIL)  simulation  facility. 


The  digital  computer  configuration  contains  an  extensive  array  of  periph¬ 
eral  devices  including  disk  files,  a  tape  drive,  and  a  graphics  display 
terminal.  A  set  of  16  priority  external  interrupts  is  included  in  the 
system;  these  may  be  triggered  by  events  which  occur  on  the  analog  com¬ 
puters  or  internally  by  the  digital  program.  The  analog  computers  are 
100-V  solid-state  machines  and  both  contain  a  wide  range  of  linear,  non¬ 
linear,  and  logic  elements. 

The  laser  target  simulator  consists  of  a  neodymium  YAG  pulsed  laser  beam 
which  is  reflected  onto  a  translucent  screen  by  a  plane  mirror;  the  mirror 
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has  two  gimballed  rotational  degrees  of  freedom  about  perpendicular  axes 
approximately  in  the  plane  of  the  mirror.  Rotation  of  the  mirror  about 
these  axes  permits  the  laser  spot  to  be  positioned  any  point  on  the 
screen.  Mirror  rotational  positions  are  servo-controlled  with  input 
voltages  derived  from  either  the  analog  or  digital  simulators.  Laser 
spot  size  and  energy  density,  are  controllable  by _a  variable  position  lens 
placed  at  the  laser  generator  output.  During  a  simulation  run,  using 
the  real  laser  sensor  in  the  loop,  the  laser  spot  is  positioned  on  the 
screen  such  that  the  actual  line-of-sight  angles  to  the  target  as  seen  by 
the  missile  are  reproduced.  The  effect  of  range  closure  between  missile 
and  target  is  simulated  by  varying  the  laser  spot  size  and  energy  density. 

The  rotational  motion  simulator  consists  of  a  three-axis  electrohydraulic, 
gimballed  flight  table  controlled  by  angular  position  inputs.  It  is  used 
to  simulate  the  missile  angular  motion  experienced  by  the  G&C  hardware 
components.  The  flight  table  is  placed  so  that  the  laser  seeker  compo¬ 
nent  of  the  guidance  hardware,  when  mounted  on  the  flight  table,  is  able 
to  sense  the  laser  spot  through  the  translucent  screen;  the  latter  acts 
to  diffuse  the  laser  spot  in  a  manner  similar  to  that  produced  by  a  real 
target . 

For  the  purpose  of  partitioning  the  mathematical  models  into  digital  and 
analog  sections,  the  models  were  considered  in  the  following  groups: 

1)  Group  a  -  This  group  contained  the  missile  dynamics  and  was 
divided  between  the  digital  and  analog  machines  with  the  translational 
equations  and  aerodynamic  table  interpolations- on  the -digital  computer. 
This  took  advantage  of  the  wider  dynamic  range  of  the  digital  machine  and 
recognized  that  the  translational  equations  had  the  lowest  frequency 
content.  The  rotational  equations  were  assigned  to  the  analog  machine 
because  the  anticipated  rotational  frequencies  required  integration  rates 
not  attainable  in  the  time  available  to  the  digital  computer  when  operat¬ 
ing  in  real-time.  However,  the  applied  aerodynamic  moment  inputs  to  the 
rotational  equations  were  generated  on  the  digital  computer  because  of  the 
necessity  of  interpolating  in  tabulated  functions  of  two  and  three  inde¬ 
pendent  variables. 

2)  Group  b  -  This  group  contained  the  target  position  and  calcula¬ 
tion  of  missile-target  miss  distance  at  intercept  and  was  assigned 
exclusively  to  the  digital  computer  because  of  the  need  for  maximum 
accuracy  in  the  miss  distance  calculation  and  its  association  with  the 
missile  translation  equations. 

3)  Group  c  -  With  the  exception  of  the  seeker  sample-hold  detector, 
Group  c  was  assigned  to  the  analog  machines  since  these  models  contained 
the  highest  frequency  content  in  the  simulation.  The  seeker  sample-hold 
output  was  generated  on  the  digital  computer  because  it  depended  on  rela¬ 
tive  translational  displacement  of  the  missile  and  target.  Jn  addition, 
however,  the  gyro  model  was  duplicated  in  the  digital  computer  to  generate 
accurate  gimbal  angles  for  the  calculation  of  detector  sample-hold  output. 
It  was  found  that  unacceptable  errors  were  introduced  in  the  transforma¬ 
tion  of  missile-target  displacement  to  seeker  axes  if  the  seeker  angle 
were  obtained  by  analog-digital ^conversion_of_analog_outputs  because  of 
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the  small  values  over  a  large  part  of  the  flight  range.  Figs.  2  and  3 
show  the  digital  and  analog  program  in  block  diagram  form. 


Fig.  2  Block  diagram  of  the  digital  Fig.  3  Block  diagram  of  the 

section  of  the  simulation.  analog  sect  ion  of  the  simulation . 


To  insure  that  the  simulation  system  was  operating  satisfactorily,  two 
dynamic  tests  were  performed.  The  analog  results  from  these  tests  were 
compared  with  results  taken  from  an  all-digital  simulation  and  a  Continu¬ 
ous  System  Simulation  Language  (CSSL)  modular  digital  simulation.  The 
procedure  for  conducting  the  test  and  interpreting  the  results  is  shown 
as  follows.  Test  No.  1  provides  a  dynamic  check  on  the  following  varia¬ 
bles:  6p  (roll  vane),  6Q  (pitch  vane),  6R  (yaw  vane),  0s  (pitch  gyro), 

4>s  (yaw  gyro),  0R  (pitch  roll  rate  gyro),  and  i</R  (yaw  roll  rate  gyro). 

The  procedure  is  to:  (1)  place  an  IC  on  the  following  variables:  P,  Q, 

R,  Wy,  and  Wz;  (2)  set  DP  *  DQ  *  DR  *  0;  and  (3)  set  control  lines  via 
the  Digital-to-Analog  Converter  (DAC) .  Test  No.  2  provides  a  dynamic 
check  on  the  following  variables:  P  (roll  body  rate),  Q  (pitch  body 
rate),  R  (yaw  body  rate),  <f>  (roll  Euler  angle),  0  (pitch  Euler  angle), 

(yaw  Euler  angle) ,  a  (pitch  angle  of  attack) ,  and  3  (yaw  angle  of 
attack).  This  was  done  in  the  same  manner  as  Test  No,  1.  Fig.  4  presents 
highlights  of  the  results  of  the  method.  For  flight  test,  the  model, 
scenario,  and  flight  parameters  require  further  discussion. 
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Fig.  4  Digital  and  analog  test  results  (Test  No.  1) . 


WEAPON  SYSTEM  CONCEPT 


Army  artillery  has  been  used  as  an  effective  weapon  for  saturation  fire 
upon  area  targets  in  the  past;  however,  the  artillery  barrage  has  been 
lacking  in  capability  to  destroy  a  point  target.  ^Because  the  armored 
moving  target  was  deemed  difficult  to  kill  and  because  it  usually  required 
a  direct  hit  to  immobilize  it,  studies  concluded  that  a  new  weapon  had  to 
be  developed  to  meet  this  threat.  It  was  suggested  that  the  ultimate 
mission  of  the  new  weapon  system  would  be  to  enhance  the  indirect  fire 
capability  by  disrupting,  delaying,  disorganizing,  and  destroying  enemy 
mechanized  forces  before  and  during  the  time  they  were  engaged  by  direct 
fire  weapons.  Thus,  the  155-mro  howitzer  with  laser  semiactive  homing 
guided  munitions  was  conceived.  The  munition  itself  is  now  commonly 
called  COPPERHEAD  (CLGP)  —  Cannon  Launched  Guided  Projectile. 

The  COPPERHEAD  CLGP  design  is  predicated  on  its  intended  use  in  a  tacti¬ 
cal  situation.  The  typical  trajectory  shown  in  Fig.  5  depicts  the  launch, 
midcourse,  and  terminal  homing  phases.  The  high-g  launch  is  shown  immedi¬ 
ately  after  which  the  stabilizing  fins  are  deployed.  On  the  ballistic 
trajectory  the  mechanism  for  guidance  and  control  is  uncaged.  After 
uncaging,  target  acquisition  can  occur;  that  is,  the  laser  energy  reflected 
from  the  target  can  be  acquired  by  the  sensor  mechanism.  Once  the  target 
has  been  acquired  by  the  sensor,  guidance  can  begin.  The  guidance  problem, 
in  which  the  heading  error  is  nulled,  begins  when  acquisition  is  achieved. 
(The  heading  error  is  defined  as  the  angle  between  the  line  of  sight  to 
the  target  and  the  velocity  vector  of  the  projectile.)  When  the  heading 
error  has  been  nulled,  the  end  game  begins  wherein  the  guidance  system 
overcomes  disturbances  such  as  those  due  tovind-or  target-motion  (see 
Fig.  6  for  a  simplified  one-axis  block  diagram). 


UN CAGE  G&C 
MECHANISMS 


ACQUISITION 
OF  TARGET 
/  / 


Fig.  5  Projectile 
sequencing  and  phases. 
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Fig.  6  One-axis 
block  diagram. 


The  COPPERHEAD  projectile  (Fig.  7)  is  155  mm  in  diameter  and  4.5  ft  long, 
and  weighs  138  lb.  The  projectile  nose  is  a  blunted  conical  dome.  Aft 
of  the  dome  is  a  primary  optical  lens  bonded  to  a  detector  assembly.  The 
detector  assembly  contains  other  optics,  a  laser  detector,  and  video 
preamps.  The  detector  and  primary  lens  assembly  are  attached  to  a  coil 
form  that  contains  both  torquing  and  spin-sustain  coils  for  the  seeker 
gyro. 


OBTURATOR 


Fig.  7  COPPERHEAD 
guided  projectile. 


A  two-axis  gimballed  gyro  using  load  transfer  bearings  and  a  f,gotcha" 
sleeve,  which  engages  the  gyro  rotor  during  launch  for  launch  survivabil¬ 
ity,  are  housed  inside  the  coil  form  assembly.  These  items  carry  the 
large  loads  incurred  during  cannon  launch  and  protect  the  gimbal  bearings 
from  brinelling.  The  gyro  rotor  has  a  mirror  surface  to  reflect  laser 
energy  onto  the  detector  and  incorporates  a  magnetic  ring  that  enables 
the  gyro  to  be  precessed  or  spun  when  the  torquing  or  spin  coils  are 
energized.  Initial  spinup  for  the  seeker  gyro  is  provided  by  a  wound 
spring. 


The  projectile  electronics  package  located  behind  the  seeker  gyro  section 
consists  of  eight  printed-circuit  cards  plugged  into- a  motherboard 


47 


assembly.  The  launch  loads  induced  on  the  seeker  gyro-  section  and  the 
electronics  package  are  carried  by  the  steel  electronics  housing  that 
mechanically  interfaces  with  the  warhead  subsystem.  The  electronics 
housing  also  incorporates  a  forward  bourrelet,  which  acts  as  a  bore  rider 
during  travel  down  the  gun  tube.  An  unconventional  roll  rate  sensor  is 
mounted  on  the  forward  bulkhead  of  jthe  electronics  jpackage^  __ 

The  steel  warhead  structure  contains  a  shaped-charge  warhead,  warhead 
fuse  with  safe  and  arm  (S&A)  provisions  and  interconnecting  electrical 
cables  for  interfacing  the  guidance  subsystem  with  the  control  subsystem. 
During  early  testing,  a  telemetry  unit  having  representative  mass  and 
envelope  characteristics  of  the  warhead  is  substituted  for  the  warhead. 

The  control  housing  contains  slots  that  permit  the  fins  and  fixed  wings 
to  be  folded  inside  the  structure  during  storage  and  cannon  launch. 

Actuator  electronics,  thermal  battery,  and  a  cold-gas  storage  bottle  for 
the  controls  are  mounted  on  the  control  actuator  unit.  The  actuator 
assembly  is  held  in  the  control  housing  by  an  aft  closure,  which  also 
incorporates  the  projectile’s  obturator.  The  control  section  also  contains 
locks  which  hold  the  fins  and  wings  inside  the  structure  and  a  cam/piston 
mechanism  which  deploys  the  wings  upon  command. 

HARDWARE- IN- THE -LOOP 

Test  Setup 

The  HWIL  simulation  is  performed  on- the  MICOM  Guidance  and- Control 
Analysis  group’s  terminal  homing  guidance  evaluation  facility.  For  HWIL 
operation  it  is  necessary  for  the  seeker  to  experience  the  projectile 
angular  orientation  and  to  sense  the  target-projectile  relative  transla¬ 
tional  displacements.  This  is  achieved  by  mounting  the  hardware  seeker 
on  a  three-axis  flight  table  to  enable  the  body-inertial  Euler  angles  for 
generation  in  real-time  while  sensing  the  laser  energy  projected  onto  a 
translucent  screen  to  simulate  the  target.  The  target  spot  is  deflected 
via  mirror  and  bellow  interfaces  to  represent  the  line-of-sight  change 
between  the  target  and  projectile  and  include  the  effects  of  relative 
range  closure  between  them.  Models  of  the  airframe,  actuator,  gyro, 
seeker,  and  compensation  networks  are  required  with  the  mathematical 
generation  of  three  flight  table  angles,  two  laser  mirror  angles,  and 
aperture  settings  for  the  laser  lens  Fig.  8  highlights  the  system 
features. 

Briefly,  the  test  sequence  consists  of  securing  the  seeker  guidance 
module  to  the  three-axis  flight  table  connecting  cables,  to  the  test  moni¬ 
tor  breakout  box,  and  to  the  computer’s  interface.  Voltage  levels  and 
polarities  are  checked  at  the  initial  power-up  condition.  Simulated 
flights  for  a  prescribed  target  are  seen  against  a  variable  laser  setting 
to  simulate  range  closure  by  increasing  the  laser  energy  transmitted  to 
the  seeker  via  the  bellows  on  the  laser  simulator.  Groups  of  five  runs 
each  quantify  the  total  range  of  the  seeker  from  maximum  range  at  detec¬ 
tion  (lowest  energy  at  acquisition)  to  maximum  energy  (impact  condition). 
Finally,  orientations  of  the  seeker  in  various  roll  positions  are  made  to 
values  of  45,  90,  and  180  deg.  If  deemed  necessary  to  improve  miss 
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distance,  vane  effectiveness  is  improved  via  settings  on  the  analog 
machine,  and  g  bias  is  increased  via  an  amplifier  gain  on  the  analog. 

A  group  of  runs  in  each  of  these  conditions  may  require  60  to  100  trails. 
These  are  all  oriented  to  maximizing  performance  prior  to  flight  testing 
at  White  Sands  Missile  Range. 


Fig.  8  Hardware-in-the-loop  simulation. 


Early  Test  Results 


A  set  of  results  from  the  hybrid  simulation  is  shown  in  Figs.  9  through 
11.  These  are  given  in  the  form  of  strip  chart  outputs  showing  the 
variation  of  parameters  over  the  terminal  homing  portion  of  the  missile 
trajectory,  and  were  obtained  with  a  target  located  12,500  ft  downrange 
from  the  missile  at  the  time  of  acquisition.  Miss  distance,  normalized 
by  the  circular  error  probability,  was  0.45  ft. 

As  soon  as  the  hybrid  real-time  simulation  was  performing  satisfactorily, 
the  guidance  hardware  components  were  inserted  into  the  simulation  loop. 
Initial  experience  immediately  highlighted  several  problem  areas  which 
had  not  been  discovered  in  the  early  open-loop  hardware  tests.  Most 
prominent  of  these  included  poor  automatic  gain  switching  in  a  set  of 
log-amps  in  the  combined  seeker  detector  and  the  gyro  torquer  loop. 

A  set  of  results  from  the  HWIL  simulation,  for  the  same  target  conditions 
as  given  for  the  hybrid  simulation  results,  is  shown  in  Figs.  12  through 
14.  Normalized  miss  distance  was  0.73  ft.  Comparison  of  these  two  sets 
of  results  shows  the  inability  of  the  simplified  models  to  represent  the 
actual  hardware  in  detail  when  high  frequency  noise  is  generated.  How¬ 
ever,  the  low  frequency  fidelity  was  consistently  good. 

Optimization  of  various  hardware  parameters  with  respect  to  missile  target 
miss  distance  at  intercept  was  then  attempted.  These  parameters  included 
navigation  ratio,  gyro  rotor  speed,  range  on  entry  to  terminal  homing, 
vane  effectiveness,  and  opt ical-gimbal  coupling.  Optimization  methodology 
consisted  of  choosing  nominal  parameter  values,  performing  a  set  of  five 
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Fig*  9  Hybrid  simulation 
pitch  channel  sample- 
hold  output  (gyro  pitch 
~~  rateO*  ~ 


Fig.  10  Hybrid  simulation 
pitch  axis  Euler  angle 
0£* 


Fig.  11  Hybrid  simulation 
pitch  channel  vane  angle. 
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simulation  runs  to  give  an  average  miss  distance  and  graphing  the  miss 
distance  against  the  parameter  to  give  an  average  miss  distance  and 
graphic  the  miss  distance  against  the  parameter  to  give  the  optimum  value 
as  a  minimum  on  the  curve.  Several  optimization  results  are  given  in 
Figs.  15  and  16.  Note  that  miss  distances  are  averaged  over  five  runs 
and  expressed  in  a  form  normalized  by_the.  circular  error  probability. 

In  Figs.  15  and  16  the  parameters  not  being  optimised  were  held  at  their 
nominal  values.  Obviously,  the  overall  optimization  procedure  required 
several  iterations  because  coupling  existed  between  many  of  the  various 
parameters.  ; 


Fig.  15  Miss  distance  as  a 
function  of  gyro  rotor  speed. 


Fig.  16  Miss  distance  as  a 
function  of  optical-gimbal 
coupling  compensation. 


Optical-gimbal  coupling,  mentioned  previously,  was  found  as  a  nonlinearity 
introduced  into  the  seeker  detector  when  the  gixnbal  displacement  permitted 
the  laser  beam  to  be  received  through  the  fixed  lens  on  an  axis  offset 
from  the  center.  An  empirical  correction  was  required  for  this  in  the 
form  of  a  compensation  loop  in  the  gimbal  torquer.  Fig.  16  shows  the 
variation  of  miss  distance  with  a  resistance  value  which  determined  the 
loop  gain. 

A  total  of  over  1500  simulated  flights  were  made  on  the  two  CLGP  designs 
over  a  period  extending  from  May  1973  to  1977.  These  include  a  sample  of 
11  projectiles  by  Martin-Marietta  and  9  by  Texas  Instruments.  Many 
variations  were  made  on  the  simulation  runs  as  described  previously. 

Thus,  it  is  difficult  to  categorize  the  data  on  an  equal  basis;  however, 
certain  statistics  have  been  calculated  and  several  examples  are 
summarized  here. 
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A  total  of  213  simulated  flights  was  made  on  four  Martin-Marietta  CLGP 
rounds  designed  for  firings  at  8-km  targets  at  WSMR.  A  total  of  217 
simulated  flights  was  made  on  one  Texas  Instruments  CLGP  at  various 
target  ranges  which  eventually  was  air  released  at  Edwards  Air  Force  Base. 
A  summary  of  these  results  is  shown  in  Table  1  (normalized  by  circular 
error  probability) . _ _ _ ■ _  ______  _  _  „  _ 

Table  1 

SUMMARY  OF  CLGP  SIMULATED  FLIGHTS 

Standard 


System 

No.  of  Flights 

Mean 

Deviation 

Texas  Inst  rumen t  s :  THU 

91  at  4  km 

2.41 

2.32 

72  at  8  km 

1.82 

2.19 

54  at  12  km 

0.76 

0.68 

Total 

217 

1.81 

2.10 

Martin-Marietta:  CG-1 

36  at  8  km 

1.34 

0.69 

CG-2 

53  at  8  km 

1.36 

0.68 

CG-3 

76  at  8  km 

1.12 

0.96 

CG-4 

48  at  8  km 

2.4< 

2.21 

Total 

213 

1.6* 

1.57 

A  summary  of  all  simulated  flights  is  given  in 
circular  error  probability) . 

Table  2 

^normalized  by 

Table  2 

HWIL 

SIMULATION  RESULTS 

System  No.  of  Flights 

Standard 

Mean  Deviation 

Texas  Instruments 

650 

2.95 

2.64 

Mart in-Marietta 

808 

1.72 

1.47 

Engineering  Development  Flight 

Test  Results 

The  first  flight  demonstrations  in  engineering  development  (ED)  took 
place  in  November  1976  to  evaluate  roll  control  and  roll  loop  stability 
limitations.  Fully  guided  ED  flights  were  initiated  in  March  1977 ,  and 
a  total  of  129  rounds  had  been  flown  as  of  August  1979.  Hardware  failures 
and  incorrectly  assumed  aerodynamic  characteristics  prevented  a  success¬ 
ful  demonstration  until  the  fifth  flight.  Subsequently  the  high-g 
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capability  of  the  components  was  proven.  After  successfully  hitting  a 
7,7-km  tank  target,  a  series  of  successes  followed  at  16,  12,  11,  7.7, 
5.6,  4,  and  3  km  range.  The  accuracy  of  the  system  (Fig.  17)  has 
exceeded  the  requirement  and  all  aspects  of  the  guidance  and  control 
loop  appear  exceptional.  Eighty-six  rounds  have  achieved  direct  hits 
on  targets.  Targets  have  been  accurately  hit  while  moving  at  speeds 
up  to  20  mph  in  both  crossing  and  head-on  attitudes.  Compatibility 
with  three  different  types  of  ground  designators  was  demonstrated. 


1  —  MEDIUM  RANGE  FI  RST  SUCCESS  STATIONARY  TANK 

2  —  MEDIUM  RANGE  20  MFH  CROSSING  TARGET 

3  —  SNORT  RANGE  ST  ATI  ONAR  Y  TARGET 
J  4  —  SNORT  RANGE  20  MFH  TARGET 

/  •— MEDIUM  RANGE  20  MFH  HEAD-ON  TARGET 

t  -  MEDIUM  RANGC  20  MFH  TARGET  LONG  DESIGNATION  RANGE 
7  —  MEDIUM  RANGE  STATIONARY  TARGET  LONG  DESIGNATION  RANGE 
«  -  MINIMUM  RANGE  STATIONARY  TARGET  HAND-HELD  DESIGNATOR 
I  —  MEDIUM  RANGC  STATIONARY  TARGET 
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Fig.  17  Typical  COPPERHEAD  guidance  and  control  performance. 


VINDICATION  BY  FLIGHT 

Two  CLGP  AD  projectiles  dedicated  to  the  RPV  designation  demonstration^ 
were  CG-11  and  CG-12.  The  guidance  sections  from  both  units  were 
evaluated  in  the  HWIL  simulation  at  the  MIC0M  G&C  analysis  facility 
during  July  and  August  1975.  CG-12  operation  was  verified  without 

incident,  but  three  evaluations  at  MIC0M  were  necessary  before  satis¬ 
factory  performance  was  obtained  with  CG-11.  The  problem  encountered 
at  MIC0M  was  an  oscillation  in  the  line-of-sight  errors  while  in  the 
high  gain  nulling  mode.  This  prevented  the  unit  from  initiating  guidance. 
The  oscillation  appeared  only  in  a' small  region  under  certain  gimbal 
angle  ranges  and  gyro  slew  rates.  This  type  problem  had  not  been 
encountered  on  any  previous  units  during  AD.  Attempts  to  diagnose  the 
specific  cause  of  the  difficulty  were  unsuccessful;  however,  the  condi¬ 
tion  was  rectified  by  a  gain  reduction  in  the  torquer  loop  and  adjust¬ 
ment  of  the  optical  transfer  functions  by  changing  the  gyro  mirror 
position. 

The  seeker  gyro  for  CG-11  was  previously  flown  in  a  MIC0M  canister  that 
experienced  failure  of  the  soft  recovery  system  and  was  subjected  to  a 
hard  impact.  The  seeker  gyro  was  subsequently  rebuilt  from  parts  avail¬ 
able  inhouse.  The  electronic  housing  for  CG-11  was  previously  discarded 
during  fabrication  due  to  cracked  lands  on  the  motherboard.  The  CG-11 
guidance  unit  experienced  another  anomaly  after  it  was  hardened.  The 
flight  sequence  logic  used  to  delay  seeker  gyro  track  and  start  of 
guidance  ceased  to  function.  Interchanging  of  suspect  components  was 
unsuccessful  and  it  was  necessary  to  add  a  circuit  modification  to 
achieve  satisfactory  performance.  After  these  problems  were  corrected, 
both  units  were  integrated  into  all  up  rounds  and  subjected  to  limited 
environmental  tests  and  inhouse  acceptance  tests.  Both  units  were 
shipped  to  WSMR  on  23  September  1975. 
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The  RPV  test  contingent  arrived  at  WSMR  and  eight  dry  run  attempts  were 
made  between  24  and  28  September  1975  without  satisfactory  results. 
Problems  were  repeatedly  encountered  with  the  VEGA  system  and  were  solved 
only  after  replacement  of  the  RPV  transponders,  repositioning  of  ground 
antennas  and  recallbratlon  of  the  equipment.  Difficulties  were  also 
encountered  in  trimming  the  autopilot  and  one  engine  failure  occurred  on 
takeoff.  Though  lasing  was  accomplished  on  a  couple  of  passes,  designa¬ 
tion  performance  and  RPV  flight  profile  control  were  unacceptable. 

In  order  to  give  increased  realism  to  the  joint  demonstration,  the  deci¬ 
sion  was  made  to  switch  from  a  proposed  12  x  12  ft  square  white  target 
panel  to  an  MAS  tank  hulk.  Reflectivity  augmentation  was  approved  and 
the  target  tank  was  taped  with  green  missile  tape.  This  produced  a 
reflectivity  of  152  (compared  to  the  normal  5  to  102)  which  provided 
adequate  reflected  energy  for  acquisition  by  CLGP. 

On  29  September,  during  the  9th  dry  run,  adequate  flight  profile  control 
was  demonstrated  and  seven  lasing  passes  were  made  against  the  tank 
target.  Designation  tracking  was  good  except  for  two  problems.  The 
tracking  sight  would  periodically  jump  (as  much  as  a  degree)  and  the 
interface  between  the  side  and  bottom  of  the  optical  dome  would  pass 
through  the  laser  beam  due  to  roll  motion  of  the  RPV.  The  sight  jump 
was  demonstrated  to  be  caused  by  interference  between  VEGA  and  FPS-16 
radars  which  were  skin  tracking  the  RPV.  The  RPV  pilot  was  not  aware  of 
the  scattering  effects  produced  by  the  dome  edge  passing  through  the 
laser  beam  and  he  was  instructed  to  maintain  the  roll  attitude  during 
the  last  20  seconds  of  flight.  The  jump  problem  was  solved  by  securing 
approval  to  drop  radar  track  on  the  RPV  (required  by  WSMR  range  safety) 
during  the  lasing  leg  of  the  racetrack.  Though  an  additional  dry  run 
appeared  desirable,  the  decision  was  made  to  attempt  a  firing  on 
3  October  provided  the  above  fixes  could  be  demonstrated  in  preliminary 
passes  against  the  target.  This  approach  minimized  the  risks  placed  on 
the  RPV  by  an  additional  takeoff  and  landing  required  by  a  dry  run 
flight. 

After  the  decision  was  made  to  proceed  with  Che  RPV  demonstration, 

MGP-11  (CG-12)  was  hardened  for  flight.  Final  assembly  and  preparation 
included  the  epoxy  hardening  of  the  aft  connector  between  the  control 
and  telemetry  sections,  final  mating  of  control  section  to  the  telemetry 
section,  and  installation  of  the  aft  closure  and  obturator.  During 
final  assembly  of  MGP-11,  the  epoxy  used  to  harden  the  connector  did  not 
set  up  properly  and  it  was  necessary  to  disassemble  the  round  to  reharden 
the  connector.  During  disassembly,  several  of  the  control  section 
connector  wires  were  broken  and  could  not  be  repaired  in  the  field.  It 
was  necessary  to  replace  the  CG-12  control  section  with  the  section  from 
CG-11,  and  repeat  the  projectile  Acceptance  Test  Procedure. 

On  3  October  1975,  the  RPV  took  off  at  11:58  a.m.  and  made  a  total  of, 
nine  passes  against  the  target.  The  periodic  jumps  observed  in  the 
tracking  sight  were  eliminated  by  shutting  down  the  missile  and  RPV 
cracking  radars  but  flight  control  adequate  to  proceed  with  a  hot  pass 
was  not  obtained.  The  RPV  was  landed  at  14:27,  refueled,  and  the  VEGA 
transponder  was  replaced.  The  second  takeoff  occurred  at  15:29,  and 
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the  RFV  made  five  passes.  Designation  and  tracking  cm  the  two  lasing 
runs  made  prior  to  the  hot  pass  were  acceptable  *  flight  profile  control 
was  much  improved  9  and  no  tracking  anomalies  were  encountered  with  both 
tracking  radars  shut  off.  The  decision  was  then  made  to  go  for  a  hot 
pass. 

CLGP  was  subsequently  launched  from  the  M109A1  at  1730  MDT  on  3  October 
and  approximately  30  seconds  later  it  achieved  a  direct  hit  on  the  tank 
being  designated  by  the  RPV.  Telemetry  data  did  not  reveal  any  project 
tile  flight  anomalies.  Examination  of  the  laser  and  designator  video 
recordings  did  show  that  designation  on  the  hot  pass  was  not  as  stable 
as  the  two  preceding  dry  run  passes.  However,  after  CLGP  acquired  at 
15  seconds  into  the  flight,  it  detected  all  pulses  except  for  three 
instances,  before  impact.  CLGP  hit  on  top  of  the  turret  near  the  base 
of  the  gun  tube  and  was  pulling  down  at  impact,  which  would  produce  a 
high  obliquity  angle  at  impact.  At  this  angle  and  location,  a  high 
probability  of  kill  would  have  been  obtained.  Fig.  18  shows  the  RPV 
and  CLGP  profile  position  during  the  actual  flight  test. 
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Fig.  18  RPV  racetrack  and  profile  position  on  launch  flight 


CONCLUSIONS 


A  comprehensive  HWIL  simulation  was  developed  to  perform  design  verifi¬ 
cation  and  extensive  statistical,  stability,  and  miss  distance  performance 
analysis  not  available  from  ah  all-digital  simulation.  The  results  to 
date  are  indicative  of  the  performance  expected  from  CLGP. 

The  hybrid  simulation  development  is  time  consuming  and  costly;  however, 
it  is  an  immensely  powerful  evaluation  and  synthesis  tool.  In  addition 
to  providing  insight  into  the  weapon  system,  it  provides  a  real-world 
method  for  optimizing  guidance  parameters  virtually  on  line.  It  repre¬ 
sents  the  highest  leveling  dynamic  validation  programmable  for  a  system 
short  of  performing  an  actual  flight  test. 
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